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Surface modelling 


3.1 INTRODUCTION 


Fraditionally, the map has been regarded as an invaluable source of 
information to both the public and industry in general. Historically, the 
cartographic map sheet was, and still is, a work of art requiring input from 
modern surveving and stereo instruments. Maps can generally be divided 
Into two groups: (1) planimetric maps, which in plan represent such natural 
and artificial features as streams, lakes, site boundaries and public works; 
and (2) topographic maps, which may include some of the above, but 
essentially represent the relief or contours of the ground. 

This chapter deals initially with topographic mapping which is the basis of 
all other forms of mapping and is not necessarily confined to land surveying. 
rhe latter part of this chapter covers surface modelling, often referred to as 
Digital Ground Modelling (DGM), Digital Terrain Modelling (DTM) or 
Digital Elevation Modelling (DEM). 

Other users, such as market researchers. engineers and scientists often 
equire the two-dimensional representation of three-dimensional data. Ina 
1ormal three-dimensional system. using conventional survey methods, the 
ligital information is often lost because the map sheet only gives a visual 
oicture. If the recovery of intormation is required for position and distances 
r areas, surveys often have to be repeated. However, if an organized 
-omputer database 1s established with a three-dimensional data grid which 
lefines a set of points in three-dimensional (X. Y. Z) space. then this 
latabase can be altered and viewed at will. These data points may lie in a 
regular or irregular (random) grid where the grid data values define the 
height’ (Z-coordinate) of each point above or below the X—Y plane. Such 
‘opographic or contour mapping programs. for the drawing of isolines along 
vhich there is a constant value. have manv users as shown in Table 3.1. 

Applications include contours in topographic mapping (equal elevations), 
sohyetals for rainfall analysis, isobars for pressure distribution, isophotes 
or light intensity, and many others. 

In addition to producing conventional two-dimensional contour maps, 
here is often a wish to display isometric projections (sometimes called mesh 
r fishnet diagrams) of the spatial data, preferably with hidden line removal 
) present a smooth uncluttered view. 
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Table 3.1 Contour applications 





User Uses 
Architects Contour plans for land developments, quarnes, reservoirs, 
etc. 
Engineers (Civil, Ground contouring, digital ground modelling, stress and 
Structural, Marine, noise contours 


Environmental, etc.) 


Geographers Topographic mapping, contouring, population 
distribution, etc. 

Geologists Interpolating and contouring the earth’s surface and 
subterranean formations 

Geophvysicists Interpolating and contouring the earth’s gravitational and 
magnetic fields, mine modelling etc. 

Hydrologists Run-off and rainfall distribution 

Meteorologists Temperature and barometric pressure studies 

Scientists Temperature and pressure studies, statistical data analysis 

Surveyors Contouring from tacheometric or aerial surveys, digital 


ground models, topographic mapping, etc. 





The programs described in this chapter allow three-dimensional data, 
either on a regular grid (PC-Grid) or at irregular or random spacing 
(PC-Random) to be entered and stored on disc for subsequent plotting. 
either- as contours (PC-Contour) or as an isometric projection (PC- 
Surfplot). In all cases the user has a choice of plotting on the screen or 
directly to a plotter. In fact it is often best to preview the plot on the graphics 
screen before dumping the graphics image on the printer or plotter to ensure 
the correct orientation and magnification. Two demonstration data files are 
also provided (Glapwell and Glenmore) to check the user's output with that 
shown in the illustrations. 


3.2 PC-GRID 


This program 1s for the creation of a data file from regularly spaced data, 1.e. 
data at the corners of a square or rectangular grid, which can then be 
accessed by either PC-Contour to draw contours or by PC-Surfplot to 
display an isometric projection on screen or plotter. 

The initial computer listing in this section, PC-Grid/Data is a blank data 
file for entering vour own data and should be saved under the title 
“GRIDRDAT”. The second computer listing is PC-Grid/Demo and includes 


data for a spoi-heap at Glapwell, and should be saved under the title 
“GRIDDEMO" as it can be used as a demonstration file to check your 
program and the output to screen or plotter. 


3.2.1 PC-Grid/Data — subroutine index 





Line numbers Function 
(it) 10-80 Initialization and control 
(b) = 100-560 Selection of application 
(c) 1000-1410 Reading and storing of data 
(d) 2000-2140 Presentation of data file information 
(¢) 3000-8990 Data file space 
({) 9000-9050 Routine to clear screen 


(2) 9500-9540 Termination of program 


(a) Initialization and control Line numbers 10—80 


Phe program ts initialized and the installation file checked for the availability 
of graphics screens and colour for subsequent display. 


(hb) Selection of application Line numbers 100-560 


[he user ts given a choice of three applications: 


(1) Read DATA and present details 
(2) Read DATA and store on disc 
(3) Prepare DATA file. 


hus this program saved under the title *GRIDRDAT” can be used initially 
O prepare a data hile. then to read the data to check it is spatially correct and 
inally to store the data on disc for access by PC-Contour and PC-Surfplot. 

A note on data file preparation is included on the screen to save the user 
constantly referring back to the book, and the user is reminded that on 
completion of entering the data. the file should be saved under another name 
‘or reference and subsequent recall. Some programs require data to be 
entered from screen prompts, but should an error occur. the data often has 
to be re-entered. To overcome this problem, this program uses DATA 
statements so that the data file can be built up in several kevboard sessions if 
eCuired. 


c) Reading and storing of data Line numbers 1000-1410 


ihe data file is then read and, depending on the number of rows (NR) and 
lumber of columns (NC) in the grid. an array Z(NR,NC) is then dimensioned 
o store the data values. These are subsequently stored on disc for quick 
cess by PC-Contour and PC-Surtplot. Before running this program. the 


user Should ensure that there is space on disc to store the data file, and that 
the disc is not ‘write-protected'! 


(d) Presentation of data file information Line numbers 2000-2140 


This routine is provided as a check to allow the user to see if the data 
information has been entered correctly before saving on disc. 


(e) Data file space Line numbers 3000-8990 


Line numbers from 3000 to 8990 are available for the insertion of data as 
suggested in the data file preparation notes at lines 500-560. 


(f) Routine to clear screen Line numbers 9000—9050 


This routine is used to clear the screen and provide lines of text at the top of 
the screen giving the title of the software program and the date. 


(g) Termination of program Line numbers 9500-9540 


This is the final display at the completion of the program, and the user is 
given an option to run PC-Contour or return to BASIC command level. 


3.2.2 PC-Grid/Data — numeric variables 


DX = X-distance of grid 

DY = Y-distance of grid 

GRID% = Application selection 

I = Integer variable used in loops 
J = Integer variable used in loops 
MAX = Maximum Z-value 

MIN | = Minimum Z-value 

NC = Number of columns in grid 
NR = Number of rows in grid 
SCD“ = Monitor display. monochrome/colour. 
SCR % = Screen resolution 

Z(NR.NC) = Arrav to hold Z-values 


3.2.3 PG-Grid/Data — string variables 


DATES = Variable used to return the date of the computer's calendar/clock in 
US format MM-DD-YYYY 

FS File name for data 

Q$ Question response (YN) 

T$ = ‘*PC-Grid/Data”’ 


3.2.4 PC-Grid/Data — BASIC program 


18 REM <GRIDRDAT> Program for IBM PC and Compatibles 
12 REM Version PC-1.080, 91-01-87 

14 REM (C) Copyright 1987 P.H. Milne 

16 REM All rights reserved 

2@ REM Creates Data File for Access by PC-Contour 
25 REM Enter Name of Data File Here for Reference 
38 CLEAR 

35 PCFS="PCSCREEN” 

40 OPEN "I",#3,PCFS$ 

45 INPUT#3,SCR& 

5@ INPUT#3,SCD% 

55 CLOSE#3 

68 ON (SCR&+1) GOTO 62,64,66 

62 SCREEN @:WIDTH 48:GOTO 8@ 

64 SCREEN 1:GOTO 68 

66 SCREEN 1:GOTO 8@ 

68 IF SCD&%=1 THEN 8@ 

78 COLOR 8,8 

88 KEY OFF 

188 CLS 

110 GOSUB 98086 

120 PRINT 

138 PRINT ® THIS PROGRAM READS INFORMATION FROM * 
148 PRINT ® REGULAR SQUARE GRID POINTS INTO A " 
15@ PRINT ® DATA FILE FOR PLOTTING BY PC-Contour” 
160 PRINT 

178 PRINT " Select application required :-" 

188 PRINT 

19@ PRINT " <1> Read DATA and present details” 

200 PRINT 

218 PRINT “™ <2> Read DATA and store on disc” 

228 PRINT 

230 PRINT " <3> Prepare DATA file” 

240 PRINT 

250 PRINT " <4> Quit "; 

260 INPUT GRID&:IF GRID%<1l OR GRID%>4 THEN 1800 

270 IF GRID%&>2 THEN 299 

288 PRINT: PRINT: PRINT “ One moment ... reading DATA file " 
298 ON GRID& GOTO 10308,1038,308,9508 

388 CLS 

318 GOSUB 90820 

328 PRINT " To check on DATA file preparation" 

338 PRINT " list lines 580-560 on screen after” 
348 PRINT ™ the Ok prompt." 

358 PRINT 

378 PRINT " After entering your own data, either” 
38@ PRINT " RUN or SAVE your program immediately.” 
3980 PRINT: PRINT 

408 END 

500 REM ***** Data File Preparation Notes ***** 
518 REM Line 3088 : Data File Name 

528 REM Line 3018 : No. of Rows, No. of Columns for grid 
530 REM Line 3028 : x-distance, Y-distance 

548 REM Line 38038 : Maximum reading, Minimum reading 
550 REM LINE 3804@- onwards, Z-data values 

568 REM To generate new data file delete Lines 3000-89990 
995 REM ***** Data File Read Routine ***** 

1000 CLS:GOSUB 9806 

1028 INPUT ® Press <Enter> when Data File Ready"”;R$ 
183@ READ FS 

1640 READ NR,NC 

18@5@ READ DxX,DY 


1860 
1070 
1080 
1698 
1100 
1116 
1120 
1130 
1142 
1190 
1200 
1216 
1220 
1238 
1248 
1250 
1266 
1276 
1280 
1298 
1308 
1310 
1320 
13306 
1340 
1356 
1366 
1370 
1386 
1390 
1400 
1416 
1995 
2000 
2016 
2020 
2030 
2048 
2050 
20690 
2870 
20680 
2090 
2180 
2110 
T QS 
2120 
2138 
2140 
2995 
30008 
8995 
9808 
9810 
9820 
9830 
9840 
9050 
9495 
95080 
9518 
95208 
953¢@ 
95486 


READ MAX,MIN 

DIM Z( (NR -—- 1),(NC —- 1)) 

FOR I = @ TO (NR - 1) 

FOR J = @ TO (NC — 1) 

READ Z(I,J) 

NEXT J 

NEXT I 

NR = INT (NR) 

NC = INT (NC) 

PRINT 

PRINT " PC Data File is now being generated...” 

OPEN "0",#1,FS$ 

PRINT#1, NR 

PRINT#1, NC 

PRINT#1, MAX 

PRINT#1, MIN 

PRINT#1, Dx 

PRINT#1, DY 

FOR I = @ TO (NR — 1) 

FOR J = @ TO (NC —- 1) 

H = Z(I,J) 

PRINT#1, H 

NEXT J 

NEXT I 

CLOSE#l 

PRINT 
BEEP:PRINT “™ PC Data File ";FS$;" stored on disc." 

PRINT 

PRINT " Do you wish to check Data File (y/n)";: INPUT Q$ 
IF Q$="n" THEN 95080 

IF Q$<>"y”" THEN 1370 
GOTO 2880 
REM ***** Data File Information ***** 
CLS 
GOSUB 9802 

PRINT “® THE DATA FILE ";FS;" HAS BEEN READ" 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT “ Do you wish to run PC-Contour now":PRINT " (y/n) 


No. Rows = : 
Minimum value 
Maximum value 
Range 

1 Map row 

1 Map column 


pe] 
3 


No. Cols.= ";NC 
";MIN 

"»MAX 

"; (MAX-MIN) 
";DY/(NR-1);" Units” 
*sDX/ (NC=1) 4". Units” 


vnun t uw 2 


IF Q$="n" THEN 9582 

IF Q$<>"y" THEN 21800 

LOAD “CONTOUR",R 

REM ***** Data Statements start here ***** 
REM Lines 3008 onwards are for your File Data 
REM ***** pC-Grid Heading ***** 
TS="PC-Grid/Data® 

LOCATE 1,2:PRINT TS 

LOCATE 1,30:PRINT DATES 

PRINT ®--]- eee ee eee esses eee eee ee eee esses eeee2e2====" 
PRINT 

RETURN 

REM ***** Program End ***** 

CLs 

GOSUB 9800 

PRINT: PRINT “ Program Complete" 

PRINT: PRINT 

END 


";: INPU 
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3.2.5 PC-Grid/Demo 


To be able to run PC-Grid/Data it is necessary to prepare a data file. It is 
suggested that once the following DATA statements have been entered into 
the file PC-Gnd/Data, it is resaved under the title “GRIDDEMO” with 
changed remarks in lines 10 and 25. The data file presented is for a spoil-heap 
at Glapwell, and the resulting data file will be automatically saved on disc as 
“GLAPWELL” for subsequent access by PC-Contour and PC-Surfplot. 

The initial data statements give the file name “GLAPWELL” to be stored 
in F$, and the number of rows (NR) as 31 and number of columns (NC) as 41. 
The Total X-distance (DX) is 400 metres and the Y-distance (DY) is 300 
metres, that is Z-values are given on a 10 metre square grid. The maximum 
(MAX) and minimum (MIN) Z-values are noted as 131.5 m and 74.1 m for 
reference and advice in later selecting the contour values required to be 
plotted. Since the zero (0, 0) of the IBM PC screen is in the top left-hand 
corner, this has been taken as the zero for data generation purposes. Thus 
the data in the top left-hand corner of the plotted area is stored in Z(9,8) and 
the data for the lower right-hand corner stored in Z((NR-1),(NC-1)). For 
maximum screen usage NR should be 75% of NC, to give a 1:1 natural grid 
plot on screen or plotter. 


3.2.6 PC-Grid/Demo — data 


3088 DATA “GLAPWELL" 

3018 DATA 31,41 

380208 DATA 480,306 

3803€ DATA 131.5,74.1 

3048 DATA 78.5,81,83,86.4,88.5,92.3,93.7,95.4,96.1,96.2,95.9,95.1,9 
3.2,92.8,92.6,92.7,93.2,93.7,93.2,91 

305@ DATA 87.6,87.2,98.3,92.6,94.8,97.3,99.9,182.4,184.9,107.2,108. 
3,109.8,118.9,111.8,112.3,112.3,112.2,112,111.5,111.0,112.5 

38068 DATA 78.7,81.2,84.7,87.6,90.4,93.6,95.6,97.5,98.4,98.3,98.2,97 
-7,96.8,95.4,95.2,95.4,96.8,96.8,94.9,92.8 

3078 DATA 98.1,89.5,98.9,93.7,96.5,99.2,102.2,105.3,107.8,110.1,112 
©6,114.3,115.7,117.2,117.9,118,117.9,117.2,116.8,115,113.2 

3880 DATA 78.3,81.1,84.7,87.4,98.2,93.4,95.9,98.2,99.8,108.7,100.6, 
108.1,99.2,98.7,98.5,98.5,98.7,98.5,97.3,95.1 

38098 DATA 92.3,89.7,98.8,94,97.2,100.6,1603.5,107.1,110.1,114.3,116. 
2,118,120.1,120.9,121.2,121.5,121.8,121.5,121.2,1280.1,119 

31808 DATA 77.9,80.2,83.6,86.1,89.3,92.7,96.0,98.1,100.3,101.8,102.4 
,102.2,1801.7,161,166.8,101,100.7,1808.1,99.2,97.4 

31180 DATA 93.7,91,98.2,93.8,97.4,1801.2,1805.1,108.5,112.3,117.4,128. 
4,121.8,122.4,123.1,123.7,123.9,124,123.9,123.3,122.6,121.6 

3128 DATA 77.3,79.4,82.4,84.3,87.5,92.3,95.8,98.1,100.4,182.9,103.7 
,103.6,103.5,163.1,162.8,102.7,102.6,102.4,101.2,1861 

3130 CATA 96.2,92.1,89.3,93.2,97.4,101.9,1805.7,1809.9,115.2,126.3,12 
2.3,123.7,124.9,125.6,125.8,126.5,126.8,126.5,126.2,125.3,124.3 

3140 DATA 76.5,78.7,880.2,83.1,86.2,98.6,95.7,97.9,100.4,183,103.9,1 
04.9,1805.2,105.4,105.2,105.1,184.9,1804.7,183.5,182.4 

3158 DATA 99.8,93.7,98.1,92.3,96.5,101.6,105.9,110.2,116.1,1208.5,12 
2.1,122.9,124.8,126.2,127.6,128.8,129.9,129.7,129,128.5,128 

3168 DATA 75.3,77.8,79.1,81.7,85.1,96.1,95.2,97.8,108.4,182.9,104,1 
@5.4,106.2,1806.6,106.7,106.9,166.7,186.3,165.7,184.5 
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3178 DATA 1602.4,97.3,91.4,98.6,95.8,101.4,105.8,118.1,115.3,119.8,1 
21,121.5,122.4,124.6,126.4,128.2,131.3,131.5,131.5,131,138.5 

3186 DATA 74.9,77.3,78.8,88.6,85.1,90.2,95.1,97.9,18080.3,1802.8,1684.5 
¢105.9,107,107.9,1808.2,1808.4,188.,187.9,107.2,106.1 

3196 DATA 163.9,100.5,93.2,98.8,95.1,108.1,104.2,109.4,113,116.5,11 
8.6,119.8,120.9,122.4,124.6,126.5,129,131.3,131.5,131,136.5 

3208 DATA 74.5,77,78.3,86.3,84.9,98.7,94.9,97.9,108.3,182.7,1804.5,1 
05.9,107,1808.2,109,1869.9,118.2,109.5,108.6,107.2 

3216 DATA 105.2,102.9,97.3,91.2,94.4,98.7,162.2,106.8,118.5,113.4,1 
15.3,117.2,118.7,1280.3,122.4,124.8,127.6,138.8,131,129,127 

322@ DATA 74.1,77,78.5,88.5,85.4,92.3,94.8,98,108.5,182.8,164.4,185 
-6,107.2,1608.5,118.1,111.2,111.5,111.2,1186.4,108.6 

3236 DATA 1€6.8,164.2,108.3,93.7,92.3,97.3,1860.3,103.8,187.3,1190.1, 
112.4,114,115.9,118.6,128.5,123,126.5,130.1,129.6,127.2,125 

3246 DATA 75.4,77.2,78.9,86.9,86.8,92.6,95.2,98.5,108.7,1803.1,184.9 
,105.8,1807.4,109,119.8,112.4,112.9,112.7,112.3,111.1 

3258 DATA 107.9,184.9,102.1,96.2,92.2,95.3,98,100.4,1804.6,187.2,189 
02,111.8,114.2,117.4,119,122.1,125.3,128.4,127,124,122 

3268 DATA 76.2,77.7,79.8,82.4,87.6,92.7,95.6,98.6,168.9,183.2,165.1 
,105.9,107.6,1809.5,111.3,112.9,114.1,114.6,114.8,113.3 

3278 DATA 116.2,187.2,182.9,98.4,91.5,94.2,96.8,98.8,181.2,184.3,18 
72£1,108.9,112.6,115.8,118.7,121.2,125.1,127,124.8,123,121.2 

3288 DATA 77.1,78.2,80.2,84.6,88.1,93.2,96.1,99.3,1801.4,183.6,185.3 
,106.8,1808,189.9,112,113.2,115.7,116.2,116.2,115.7 

3298 DATA 112.8,188.1,104.6,1880.7,92.6,93.2,95.2,97.6,98.8,162.1,10 
4.8,107.8,111.1,114.9,117.8,128.7,125.1,125.1,123.8,121.5,119.3 

3306 DATA 77.8,79.6,82.1,85.7,89.8,94.1,97.4,108.2,102.4,104.3,105. 
7,107.1,168.3,1180.8,112.2,114.3,116.4,117.7,117.9,116.9 

3318 DATA 114.9,109.8,105.7,102.3,94.9,93,94.6,96.9,98.2,108.3,1863. 
7,106.9,118.1,113.8,117.6,1286.3,123.8,123.8,122.6,128.2,118 

3328 DATA 78.3,80.2,83.1,87.4,92.1,95.2,98.3,101.2,183.1,184.9,106. 
2,167.6,1808.9,111.2,112.7,114.8,117.4,119.3,128.2,118.7 

3338 DATA 116,112,1807,182.8,97.1,92.9,94,96,97.9,99.8,182.7,185.3,1 
69.6,113.2,117,119.8,122.6,122.6,120.5,118.7,116.5 

3340 DATA 78.8,82.6,85.6,89.9,93,97.4,108.2,102.4,104.2,105.6,106.9 
,108.1,109.5,111.4,113.1,115.2,118.3,120.8,121.7,120.3 

3358 DATA 116.4,112.4,107.4,1803.6,98.8,92.8,93.8,95.1,97.6,99.6,162 
-3,104.8,108.2,112.3,115.8,118.6,121.5,128.8,118.3,116.4,114.3 

3366 DATA 79.9,82.7,87.9,92.6,95.9,99.8,101.8,183.3,105.1,1806.3,107 
-6,108.8,1809.9,112,114.2,116.8,119.1,121.7,123,128.9 

3378 DATA 117.2,113.2,108.2,184.3,100.4,92.7,93.8,94.7,97.2,99.1,18 
1.5,104,1807.6,111.5,115.1,117.3,119.5,118,116,112.6,169.5 

338@ DATA 81.2,85.7,98.8,94.8,97.9,100.9,102.8,104.8,165.8,107.1,18 
8.2,1809.4,118.5,112.9,115.2,117.7,119.7,122.2,123.3,121.1 

3398 DATA 117.4,113.5,109.3,105.2,101.8,93.5,94,94.6,96.5,98.5,181, 
183.2,107,110.1,113.7,115.9,117.2,115.1,112.6,118.2,108 

3408 DATA 82.4,87.6,92.6,96.1,98.7,102.1,1803.3,105.6,106.7,1868.3,10 
8.8,109.9,111.6,113.7,116.3,118.2,128.1,122.7,123.3,128.8 

3416 DATA 117.4,114.1,189.8,105.9,182.4,94.7,94.2,94.6,95.9,98,180. 
1,162.6,105.9,1809.8,112.3,113.6,113.8,112.7,118.5,189,107.5 

3420 DATA 84.9,89.5,93.2,97.3,100.4,1802.9,1804.8,1806.2,107.4,188.4,1 
@9.4,1186.6,112.4,114.3,116.9,118.8,120.6,122.9,123,1280.4 

3438 DATA 117.2,114.8,110.1,106.2,182.4,95.2,94.3,94.5,95.8,97.8,99 
-6,101.8,105.1,168,118.8,111.3,110.7,109.7,108.3,107.2,1606.1 

3448 DATA 87,90.6,94.8,98.8,101.6,103.9,106,107.3,108.1,169,169.9,1 
11.7,113,115.8,117.4,119,128.8,123.2,122.6,119.8 

345@ DATA 117.1,114.9,110.2,1807,183,96.3,94.5,94.5,95.1,96.8,97.9,9 
9.3,102.2,104.8,106.3,1806.9,106.9,106.3,186.3,1805.1,183.9 

3468 DATA 88.1,92.6,97.2,180.3,182.9,105.3,166.9,108.2,188.9,169.7, 
118.8,112.1,114.2,116.3,118.1,119.2,121.4,123.8,121.7,118.7 

3478 DATA 116.8,115.1,111.3,1607.6,183.7,99.8,95.1,94.6,94.7,95.5,96 
-8,97.4,99.3,1800.7,101.9,1802.6,182.9,183.7,183.8,183,162.2 

3486 DATA 89.9,93.5,97.6,101.2,103.7,1096.1,167.4,108.4,189.4,118.3, 
111.2,112.7,114.9,116.6,118.7,119.5,121.7,123.4,128.6,118.2 


3490 CATA 116.2,114.8,111.9,108.5,104.6,100.7,97,95.1,94.8,95.3,96. 
1,96.7,97.2,97.8,98.9,99.7,1808.1,101.9,182.8,103,103.5 

3508 DATA 92,95.1,99.2,102.4,104.8,106.6,107.9,108.8,109.7,118.6,11 
2.1,113.3,115.2,116.8,117.8,119.5,121.7,121.7,119.6,117.4 

3516 DATA 116,114.6,112.1,169.7,105.2,101.3,99.2,97.1,94.9,95.1,95. 
4,95.8,96.2,96.6,97,97.2,97.4,99.8,1801.6,162.7,183.8 

3526 DATA 93.2,96.8,160.7,103.1,1865.8,1607.2,108.2,129,118.1,111.2,1 
12.4,114.2,115.9,116.9,117.4,119.2,128.8,126.1,118.3,116.8 

3538 DATA 115.9,114.3,112.1,118.8,107.4,162,160.1,98.2,96.5,97.6,97 
02,97.3,97.4,97.2,97,96.5,96,95.5,97.8,1808.2,182.6 

3548 DATA 93.7,97.4,101.8,1804.2,106.1,167.4,1898.5,1209.5,118.5,111.6 
2112.7,114.5,116.1,117,117.3,118.5,119.2,118.5,117.4,116.3 

3558 DATA 115.6,113.8,111.8,1186.8,1808,1802.7,100.2,97.6,97.3,160.1,9 
9.8,99.7,99.5,98,97.6,97.2,96,96,97.2,98.8,188.4 

3566 DATA 93.7,97.6,1802.3,104.8,106.3,167.8,198.8,1809.6,1180.7,111.9 
0113.1,115.1,116.3,117.3,117.3,117.6,118.6,118,117,116 ‘ 

3576 DATA 114.9,112,1186.9,110.2,168,105.1,102.4,99.3,99.3,188.5,1808 
-8,100.8,1880.7,180.4,99.5,98.8,97.6,96.2,96.3,97.4,98.5 

35886 DATA 93.7,97.7,162.3,104.8,1806.6,108,108.9,169.7,111.2,112.2,1 
13.5,115.1,117.4,118.1,118,117.2,118.3,117.9,117,115.8 

3598 DATA 112.1,1809.6,1809.2,108.5,1867.9,104.9,102.7,161.8,181.8,1982 
©1,102,162,1801.9,1801.7,1808.8,99.9,98.7,97.3,96.4,96.5,96.6 

3686 DATA 92.6,97.4,162.1,104.6,186.9,1607.9,108.8,1089.8,111.4,112.3 
2113.7,116,119,1280.7,121,120.4,117.3,117.8,117,115.6 

36186 DATA 110.9,109.4,168.8,1808.2,167.3,1804,103.5,163.3,103,102.9,1 
62.8,102.7,102.6,182.2,1801.9,161.7,99.8,98.2,97.3,96.8,96.3 

36280 DATA 92.3,95.7,101.4,1804.1,166.5,187.9,198.8,109.9,111.4,112.4 
2114.1,116.8,119.8,121.5,121.8,121.4,117.9,117.9,117.4,114.8 

3638 DATA 111.2,1809.5,1808.5,197.4,105.2,1804.6,103.6,103.5,183.5,103 
-5,103.5,183.3,1803.1,1802.8,182.3,101.5,106.8,99.1,98.2,97.1,96 

3648 DATA 92,96,100.7,1603.6,106.1,107.7,108.6,109.7,111.6,112.5,114 
©5,117.5,118.9,121.5,123,122.4,118.5,118,117.8,115.2 

3658 DATA 111.5,109.6,108.3,106.6,1603.2,104.3,103.7,1803.6,1604,104.1 
2184.2,103.9,1803.6,103.4,103.7,101.3,101.3,99.7,99.1,97.4,96.2 


3.3 PC-RANDOM 


Not all spatial data is collected or sampled on a regular grid basis and this 
program allows contours to be drawn from irregularly or random spaced 
data. The data has first to be read by the program PC-Random/RData and 
then a disc data file must be prepared for access by PC-Contour or 
PC-Surfplot. 

The initial computer listing in this section, PC-Random/RData is a blank 
data file for entering your own data and should be saved under the title 
“RANDRDAT”. The second computer listing is PC-Random/Demo and 
includes data for a tacheometry survey at Glenmore in Scotland; it should 
be saved on disc under the title “RANDDEMO” as it can be used as a 
demonstration file to check your program and the output to screen or 
plotter. 


3.3.1 PC-Random/RData — subroutine index 





Line 
numbers Function 
(a) 10-90 Initialization and control 


(b) 100-560 Selection of application 

(c) 1000-1540 Reading and storing of data 

(d) 2000-2440 Presentation of data file information 
(e) 2500-2990 PC-Random iteration routine 

(f) 3000-3490 Graphics screen display 

(g) 4000-4850 Plotter routines 

(h) 6000-6270 Printout of data file details 

(i) 9000-9050 Routine to clear screen 

(j) 9500-9650 Termination of program 


(a) Initialization and control Line numbers 10-90 


The program is initialized and the installation file checked for the availability 
of graphics screens and colour for subsequent display of the random data 
points. 


(b) Selection of application Line numbers 100-560 
The user is given a choice of three applications: 


(1) Read DATA and present details 
(2) Read DATA and store on disc 
(3) Prepare DATA file 


Thus this program stored under the title “RANDRDAT” can be used 
initially to prepare a data file, then to read the data to check it is spatially 
correct, and finally to store the data on disc for access by PC-Contour or 
PC-Surfplot. 

To save the user constantly referring back to the book, a note on data file 
preparation is included on the screen, and the user reminded that on 
completion of entering the data, the file should be saved under another name 
for reference and subsequent recall. Data is entered in data statements as 
described in 3.2.1(b). 


(c) Reading and storing of data Line numbers 1000-1540 


The data file is then read, first its name, F$ and then the X-coordinates of 
the left-hand side and right-hand side of the plot, followed by the bottom 
and top coordinates, i.e. the plotting area is defined by X-min, X-max, 
Y-min, Y-max. The number of data points, N, is then read and arrays set 
up to store the X, Y and Z-coordinates of the points. There is no need to 
enter the maximum and minimum Z-values in this program as this routine 
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Fig. 3.1 PC-Random/Demo plot of data points for Glenmore tacheometry survey on 
Screen 1 (medium-resolution). 


automatically locates them and stores them in the variables MAX and MIN 
respectively. 

The plotted area is then split up into a grid; the maximum number of 
columns is taken as 41 and the maximum number of rows is taken as 31. If, 
for example, a square area was entered, the number of columns would be 
automatically reduced to 31 to ensure the correct scaling of the data on the 
screen or plotter. As PC-Contour and PC-Surfplot require grid points for the 
interpolation of contours, these grid points are calculated in the routine at 
line 2500 before being stored on disc. The user should ensure before running 
this program that there is space on disc to store the data file, and that the disc 
is not ‘write-protected’. The user should also be patient during this 
interpolation process, as with a microcomputer running at 4.77 MHz, the 
time for the data file generation can take 45-50 minutes. Faster machines 
will take less time. 


(d) Presentation of data file information Line numbers 2000-2440 


This routine is provided as a check to allow the user to see if the data 
information has been entered correctly before storing on disc. The user is 


PC-RANDOM 89 


208.3 ~ abd 
al2.2 204.1 197.4 


wc - 196.2 200.9 
193.8 
T 86,5 et 


202.4 192.9 
199.1 182.6 


196.2 187.2 
mt 193.7 91.3 178.5 185.6 


> 199.9 1036 177.2 
186,1 176.8 


188.9 
182.9 182.8 182.3 174.3 
178.2 178.3 178.5 173.3 


190.4 © 186.5 7.7 


187.3 





Fig. 3.2 PC-Random/Demo plot of ground levels for Glenmore tacheometry survey 
on Screen 2 (high-resolution). 


also given an opportunity to display the data on the graphics screen or 
plotter. 


(e) PC-Random iteration routine Line numbers 2500-2990 


As mentioned in (c) this iteration routine is used to scan the PC-Random 
data and calculate the value of each grid point from the five nearest data 
points using an inverse square distance weighted interpolation. As in the 
normal hand drawing of contours, the final contours will only be as good as 
the selection of the data points to accurately represent the ground surface. 
Students should resist the temptation to reduce data coverage on fairly level 
ground as this will impair the accuracy of the final contours. 


(f) Graphics screen display Line numbers 3000-3490 


The display screen chosen for the display of either the data points or their 
Z-value is shown in Fig. 1.4. The screen is first set up to cover X-min, X-max, 
Y-min, Y-max and then the position of each data point calculated and 
labelled in text mode, using the variables CROW and CCOL for the rows 
and columns. On completion, the user can obtain a graphics printer screen 


dump or plot. Hard copies of both the graphics screen and plotter can be 
obtained through this menu of either the data points, Fig. 3.1, or their 
Z-value, Fig. 3.2. 


(g) Plotter routines Line numbers 4000-4850 


Three plotter routines are provided for the Apple 410, Epson HI-80 and 
Hewlett-Packard HP-7475A colour plotters, as discussed in Chapter 2. The 
user has a choice of plotting either the data points or their Z-values? On 
completion the plot is labelled with the data file name and date of plotting. 


(h) Printout of data file details Line numbers 6000-6270 


This printout routine is provided to keep a record of the data information 
entered, X-min, X-max, Y-min, Y-max, the number of points and their_X, 
Y, Z-coordinates for future reference. 


(i) Routine to clear screen Line numbers 9000—9050 


This routine is provided to clear the screen before each plot and provide lines 
of text at the top of the screen to identify the program title and data file 
presentation, together with the current date. 


(j) Termination of program Line numbers 9500—9650 


This is the final display at the completion of the program, and the user is 
given an option to run PC-Contour or return to BASIC command level. 


3.3.2 PC-Random/RData — numeric variables 


CCOL = Column number for text display on graphics screen 
CROW = Row number for text display on graphics screen 
DX = X-distance of plotted area 

DY = Y-distance of plotted area 

FX = X-distance between grid points 

FY = Y-distance between gnd points 

GPX8-2 = Variables to plot grid nodes on screen 

GPY8-2 = Variables to plot grid nodes on screen 

GX = X-coordinate at lower left-hand corner of plot 
GY = Y-coordinate at lower left-hand corner of plot 
HP% = Integer variable to identify plotter 


I = Integer variable used in loops 


J = Integer variable used in loops 

K = Integer variable used in loops 

LI-LS = Variables to identify level points in grid 
MAX = Maximum Z-value in data 

MF = Ratio of DX to DY 

MG = Multiplication factor for plotted display 


x «@ 
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Multiplication factor for horizontal resolution 
Minimum Z-value in data 

Multiplication factor for vertical resolution 
Multiplication factor for Epson HI-80 plotter 
Number of data points — 

Number of columns in grid 

Number of rows in grid 

Screen display variable 

X-coordinate at upper right-hand corner of plot 
Y-coordinate at upper right-hand corner of plot 
Array for data point identification 

Integer variable for plotter selection 

Reduced level of data output 

Array to store reduced levels 

Integer variable to select application 

Integer variable to select line colour 

Variables to calculate grid node levels 

Integer variable to define screen display 
Integer variable to select screen or plotter 
Select display, identification number or level of data point 
Integer variable to select screen or plotter 
Integer variable for screen resolution 

Array to store identification number or level of data point 
Horizontal display width in user units 

Vertical display height in user units 
X-coordinate of grid point 

X-coordinate of left-hand side of plot 
X-coordinate of right-hand side of plot 

Array for X-coordinates of data points 
Y-coordinate of gnd point 

Y-coordinate at foot of plot 

Y-coordinate at top of plot 

Array for Y-coordinates of data points 
Z-coordinate of interpolated gnd points 

Array for Z-coordinates of data points 





3.3.3 PC-Random/RData — string variables 


DATE$ 


Variable used to return the date of the computer’s calendar/clock in 
US format MM-DD-YYYY. 
File name 
Entry for graphics dump or printer 
“PC-Random” 
“PCSCREEN”, installation file for graphics 
Question response, (y/n) or (Y/N) 
P$ + F$ 


F$ 
GD$ 
P$ 
PCF$ 
Q$ 
T$ 
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3.3.4 PC-Random/RData — BASIC program 


560 


REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 


CLEAR 


<RANDRDAT> Program for IBM PC and Compatibles 

Version PC-1.8, 81-01-1987 

(C) Copyright 1987 P.H. Milne 

All rights reserved 

Creates Data File for Access by PC-Contour 

Plot available of DATA on Apple, Epson or Hewlett-Packard 
Setting in Line 48098 "COM1:1200,N,7,1,CS2008,DS15800" 
Storage space for DATA in lines 588@-onwards 


PCFS=" PCSCREEN® 


OPEN 


“1",#3,PCFS 


INPUT@3,SCRE 

INPUT#3,SCD% 

CLOSE#3 

ON (SCR8+1) GOTO 62,64,66 
SCREEN @:WIDTH 4@:GOTO 88 
SCREEN 1:GOTO 68 

SCREEN 1:GOTO 88 

IF SCD%=l1 THEN 828 

COLOR 8,98 

KEY OFF 


FS=" 
CLS 


GOSUB 98600 
PRINT 
PRINT ® THIS PROGRAM READS INFORMATION FROM” 
PRINT " RANDOMLY SPACED POINTS INTO A® 
PRINT °® DATA FILE FOR PLOTTING BY PC-Contour" 
PRINT 
PRINT " Select application required :-" 
PRINT 
PRINT “" <l> Read DATA and present details" 
PRINT 
PRINT * <2> Read DATA and store on disc” 
PRINT 
PRINT " <3> Prepare DATA file" 
PRINT 
PRINT ° <4> Quit "*; 
INPUT RAND%:IF RAND%<1l OR RAND%>4 THEN 1902 
ON RAND% GOTO 10308,1830,3¢€08,9500 


CLS 


GOSUB 98086 

PRINT " To examine Data File preparation” 
PRINT " list lines 506-566 on screen after" 
PRINT " the Ok prompt" 

PRINT 

PRINT " After entering your own data, either” 
PRINT ® RUN or SAVE your program immediately.” 
PRINT: PRINT 


END 
REM 
REM 
REM 
REM 
REM 
REM 
REM 


week Data File Preparation Notes ***** 


Plotter 


Line 5868 : Data File Name 

Line 58@18 : X-coordinates LHS, RHS 

Line 5828 : Y-coordinates bottom,top 

Line 5038 : Number of data points 

Line 5048 - onwards : DATA Point No, X-coord, Y-coord, Z-coord 
To create own data file start at line 50086 


1800 CLS:GOSUB 98988 

1020 INPUT " Press <Enter> when Data File Ready °;R$ 
1030 READ FS$ 

1048 READ XL,XR 

1858 REAC YB,YT 


1060 
1870 
168¢ 
16986 
1188 
1118 
1120 
1136 
1148 
1150 
1168 
1178 
1188 
1190 
1206 
1218 
1220 
1230 
1246 
1250 
1266 
1270 
12880 
Os: J 
12908 
1300 
BEING 
1310 
1320 
1336 
1348 
1350 
1366 
1370 
1382 
1390 
1480 
1410 
1420 
1430 
145@ 
1462 
147@ 
1480 
1580 
151¢ 
1526 
153¢ 
154@ 
1995 
2000 
2262 
2004 
20806 
2008 
2018 
2820 
2830 
2048 
2850 
2868 
2870 
2089 
28990 
2180 
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READ N 
DIM P(N) ,X(N),¥(N),Z(N),R(N) ,S(N) 
PRINT: PRINT 
PRINT " DATA now being read from file"® 
MAX=-25008 
MIN=25000 an 
FOR I=1 TO N 
READ P(I) ,X(I),Y¥Y(I) ,2(1I) 


IF Z(I) > MAX THEN MAX = Z(TI) 
IF Z(I) < MIN THEN MIN = Z(TI) 
NEXT I 

DX = XR - XL 

DY = YT - YB 

MF = DX / DY 

NC = 41 


NR=(NC-1)/MF + 1 

IF NR<=31 THEN 1252 
NR = 31 
NC=(NR-1) *MF + 1 
NC = INT (NC+.1) 
NR = INT (NR+.1) 

FRINT: IF RAND%=1 THEN 1518 
PRINT " Has Data been stored on Disc already":INPUT " 
F QS="y" THEN 15186 

IF Q$<>"n" THEN 1278 
PRINT: PRINT:PRINT “ PC-Contour Data File :";F$:PRINT ” 
GENERATED .....-” 

OPEN “0O",#1,FS$ 

PRINT#1, NR 

PRINT#1, NC 

PRINT#1, MAX 

PRINT#1, MIN 

PRINT#1, DX 

PRINT#1, DY 
FX = (DX / (NC = 1)) 
FY = (DY / (NR - 1)) 


FOR I = @ TO (NR - 1) 
Y= YT- I * FY 

FOR J = @ TO (NC - 1) 
X = XL + J * FX 

GOSUB 2508 
°*NEXT J 

NEXT I 

CLOSE#1 


BEEP:PRINT:PRINT “ Data File -";FS$;" - Closed” 
BEEP:PRINT:PRINT “® Do you wish Printout of Data (y/n)"; 
INPUT QS 

IF QS="y" THEN 6000 

IF Q$<>"n" THEN 1518 

REM *****® Data Information and Plot Selection ***** 
CLS 

SCREEN SCR&:IF SCR¥=8 THEN 20910 

IF SCR¥=1 THEN 2808 

SCREEN 1:GOTO 28018 

COLOR 8,0 

GOSUB 9688 

PRINT " THE DATA FILE “;F$;" HAS BEEN READ" 

PRINT 


PRINT " No. Rows = ";NR; ® No. Columns = *°;NC 
PRINT " Minimum Value = "°;MIN 

PRINT * Maximum Value = ";MAX 

PRINT * Range = "; (MAX-MIN) 

PRINT " 1 Map Row =z "sDY/(NR-1);" Units" 
PRINT " 1 Map Column = *;DX/(NC-1);" Units® 


BEEP 


(y/n) 


IS NOW 
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2110 PRINT: PRINT " Select Display :-" 

2128 PRINT: PRINT " <1> SCREEN, <2> PLOTTER or <3> QUIT": 
2138 IF SCR&<>8 THEN 2158 

21480 PRINT: PRINT "° N.B. Option <l> not available” 

215@ INPUT SEL% 

2166 IF SEL&<l OR SEL%>3 THEN 2188 

2178 IF SCR%=@ THEN IF SEL&=] THEN 2008 

2188 ON SEL% GOTO 2308,48008,9508 

2295 REM Choose Screen for Plotting 

2308 CLS 

2382 SCREEN SCR 

2304 IF SCR%=1 THEN 2308 

2306 SCREEN 1:GOTO 23106 

2388 COLOR 8,9 

231@ GOSUB 9800 

2320 PRINT:PRINT TAB(3);"Select SCREEN";TAB(2@);"<l> Colour -" 


2330 PRINT TAB(2@) ;* Low resolution" 
2340 PRINT: PRINT TAB(20);"<2> Black & White -” 
2358 PRINT TAB(26);” High Resolution" 


2360 PRINT: PRINT TAB(28);: INPUT SR% 

2370 IF SR&¥<1 OR SR¥>2 THEN 23298 

2380 PRINT:PRINT:PRINT TAB(3);"Select Data for Screen :-" 
2398 PRINT:PRINT TAB(20);"<1> Point No. or” 

2400 PRINT:PRINT TAB(20);"<2> Z-Value ";: INPUT SD% 
2410 IF SD&<1 OR SD%>2 THEN 2388 

24280 IF SR*=1 THEN 30880 

2430 SCREEN 2 

2448 GOTO 3800 

2495 REM ***** pC-Random Iteration Subroutine ***** 
2508 FOR K = 1 TON 

251@ R(K) = (X - X(K)) ~ 2 + (¥Y - Y(K)) ~*~ 2 

2520 IF R(K) = @ GOTO 298@ 

253@ NEXT K 

2548 RL = 25000 

2558 FOR K = 1 TON 

256@ IF R(K) > RL GOTO 2598 

2570 Ll = K 

2580 RL = R(K) 

2598 NEXT K 

2608 RL = 25002 

2610 FOR K = 1 TON 

2620 IF K = Ll GOTO 26628 

2630 IF R(K) > RL.GOTO 2662 


2648 L2 = K 

265@ RL = R(K) 

2660 NEXT K 

2678 RL = 25808 

2680 FOR K = 1 TON 

26980 IF K = L1 GOTO 2748 
2700 IF K = L2 GOTO 2748 
2718 IF R(K) > RL GOTO 2748 
2720 L3 = K 

2730 RL = R(K) 

2740 NEXT K 

2758 RL = 2500@ 

2760 FOR K = 1 TON 

2778 IF K = Ll GOTO 2838 
2786 IF K = L2 GOTO 2838 
2798 IF K = L3 GOTO 2838 


2800 IF R(K) > RL GOTO 2838 
2810 L4 = K 

2820 RL = R(K) 

2830 NEXT K 

2840 RL = 25008 

2858 FOR K = 1 TON 


2868 
2870 
2880 
2896 
2960 
29186 
2928 
2938 
2940 
4) + 
2950 
2960 
2978 
2988 
2996 
2995 
3800 
30106 
38208 
3630 
3848 
3058 
3060 
30708 
3080 
3898 
3100 
31190 
3128 
3130 
3146 
3158 
3160 
3176 
3188 
319@ 
3195 
3200 
3218 
3215 
3220 
3236 
3246 
3250 
3255 
3256 
3260 
3265 
3278 
32890 
3290 
3295 
3300 
33186 
3320 
3330 
3348 
3350 
3360 
3378 
3480 
341¢ 
3428 
343¢ 
3448 
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‘Mm 


IF K = Ll GOTO 29380 
IF K = L2 GOTO 29380 
IF K = L3 GOTO 2938 
IF K = L4 GOTO 2938 
IF R(K) > RL GOTO 2938 
LS = K —_— 
RL = R(K) 
NEXT K 
Sl = Z(L1) / R(L1) + Z(L2) / R(L2) + 2(L3) / R(L3) + Z(L4) / R(L 


Z(L5) / R(LS) 

S2 = 1/ R(L1) + 1 f/ R(L2) + 1 / R(L3) + 1 S R(L4) + 1 S R(LS) 
ZP = Sl / S2 

GOTO 2998 

ZP = Z(K) 

PRINT#1, ZP:RETURN 
REM ***** Screen Plot *#**## 
CLS 
GX=5*SR&% sGY=5 : MG=6 : MH=0 :MV=08 :SX=240*SR%:SY=180: PP=90 
MH=SX/ (DX*SR&) :MV=SY/DY 

IF MH<=MV THEN MG=MH ELSE MG=MV 
WINDOW (80,8) ~ (320*SR%, 200) 

IF SR%=1 THEN RGB%=2 ELSE RGB%=3 
PX=GX+MG * DX*SR¥% : PY=GY+MG * DY 

LINE (GX,GY)~(PX,PY) ,RGB&,B 

IF SR¥=]1 THEN 31060 

LINE (GX-1,GY)~-(PX-1,PY) ,RGB&,B 
FOR I=10 TO (NC-10) STEP 16 
FOR J=10 TO (NR-16) STEP 160 
GPX0=GX+(6*I*SR%) 
GPX1=GPX8~5*SR& :GPX 2=GPX8+5*SRE 
GPY@=GY+(6*J) 
GPY1=GPY@0-S5:GPY2=GPY6+5 
LINE (GPX1,GPY8@) =~ (GPX2,GPY6) ,RGB% 
LINE (GPX8,GPY1) —- (GPX8@,GPY2) ,RGB& 
NEXT J 
NEXT I 

IF PP=1 THEN 33080 

* Start of Screen Plot 

LOCATE 1,2:PRINT “PC-Random/RData ~ ";FS$ 
LOCATE 1,30*SR&:PRINT DATES 
FOR J=1 TO N 
IF SD%=1 THEN S(J)=P(J) 

IF SD%=2 THEN S(J)=Z(J) 
CCOL=(GX+((X(J)—-XL) *MG) *SR&) /8: CROW= (GY+(Y(J)-YB) *MG) /8 
IF CCOL<2 THEN CCOL=2 

IF CCOL>28*SR%& THEN CCOL=28*SR% 

IF CROW<=3 THEN CROW=3 

IF CROW>22 THEN CROW=22 

LOCATE (26-CROW) ,CCOL-1 
PRINT S(J) 
NEXT J 
PP=1:GOTO 30786 
BEEP 
LOCATE 3,32*SR&:PRINT "Screen" 
LOCATE 4,32*SR&:PRINT “Dump(y/n) ® 
LOCATE 5,32*SR%&:PRINT " ™ 
LOCATE 5,32*SR&: INPUT GDS 

IF GDS="n" THEN 3408 

IF GD$<>"y" THEN 3360 
LCOPY 
BEEP 
LOCATE 9,32*SR&:PRINT “Select” 
LOCATE 18,32*SR&:PRINT "Next” 
LOCATE 11,32*SR&:PRINT "Choice ;* 
LOCATE 13,32*SR&:PRINT °<1l>Screen” 
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3458 LOCATE 14,32*SR&:PRINT "<2>Plot” 

3468 LOCATE 15,32*SR&:PRINT "<3>Quit” 

3465 LOCATE 16,32*SRt:PRINT ” 

3478 LOCATE 16,32*SR&: INPUT SELECTION 

3480 IF SELECTION <1 OR SELECTION>3 THEN 3400 

3498 ON SELECTION GOTO 2308,4800,9508 

3995 REM *****® Colour Plotter Display ***** 

4008 GX=8:GY=6:SX2192:SY2#1 20 

4018 MG=@8:MH=6 :MV2@6 

4028 MH=SX/DX:MV=SY/DY:MX=2678/SX 

4030 IF MH<=MV THEN MG=MH ELSE MG=MV 

4040 PX=MG*DX:PY#=MG*DY 

4858 CLS:GOSUB 998080 

4051 PRINT:PRINT * Select Plotter"™;TAB(18);"<1> Apple 416" 

4052 PRINT:PRINT TAB(18);"<2> Epson HI-80":PRINT:PRINT TAB(18);"<3> H 
ewlett—Packard";:INPUT PLOT’ 

4054 IF PLOT&<1 OR PLOT%>3 THEN 40658 

4056 IF PLOT%=1 THEN HP=3:REM Apple 

4058 IF PLOT&=2 THEN HP=4:REM Epson 

4059 IF PLOT%=3 THEN HP=1]1:REM HP-GL 

4068 PRINT:PRINT:PRINT TAB(3);"Select Data for Plotter :-" 

4062 PRINT:PRINT TAB(18);"<1> Point No. or” 

4064 PRINT:PRINT TAB(18) ;"<2> Z-Value "; 

4066 INPUT SD&% 

4068 IF SD&<1 OR SD&>2 THEN 48680 

4078 PRINT:PRINT: PRINT “ Place paper on Plotter and switch on" 
4080 PRINT: INPUT ° To Continue Press <Enter> ";RS$ 

4098 OPEN "COM1:1200,N,7,1,CS28008,DS158000" FOR OUTPUT AS #2 

4100 ON HP GOTO 4250,42508,4185,4208 

4103 REM ***** Apple 410 ****#*# 

4105 PRINT#2,"SP3;”" 

4110 PRINT#@2,"WD -—5,-5,192,125;" 

4120 PRINT#2,"PS1;" 

4140 PRINT#2,"MA";GX;",";GY;”" 
414) PRINT#2,"DA";PX;",";GY;" 
4142 PRINT#2,"DA";PX;",";PY;" 
4143 PRINT#2,"DA";GxX;",";PY;" 
4144 PRINT#2,"DA";GxX;",";GY;";3" 

4150 FOR I=18 TO (NC-180) STEP 18 

4155 FOR J=198 TO (NR-1@) STEP 104 

4160 PGXO0=GX+I*PX/ (NC—1) 

4161 PGX1=PGX0~2: PGX2=PGX6+2 

4162 PGYO=GY+J* PY/ (NR-1) 

4163 PGY1l=PGY@-—2: PGY2=PGY0+2 

4170 PRINT#2,"MA";PGX1;",";PGY0;";":sPRINT#2, "DA"; PGX2;",";PGY8;";" 
4175 PRINT#2,"MA";PGXO;",";PGY1;";"sPRINT#2,"DA" ;PGXO;",";PGY2;";" 
4180 NEXT J 

4185 NEXT I 

4198 GOTO 4308 

4195 REM ***** Epson HI-80 ***** 

4200 PRINT#2,"IW0,8,2678,1920" 

4201 PRINT#2,"MA";GX;",";GY 

4202 PRINT#2,"0OR” 

4283 PRINT#2,"SP1":PRINT#2, "LTO" 

4204 PRINT#2,"MA";:GX;",";GY 

4205 PRINT#2,"DA"sPX*MX;3",";GY;", "3 PX*MX3","sPY*MX;","3GX3",";PY*MX;" 
0" GGX;",";GY 

421@ FOR I=18 TO (NC~10) STEP 18 

4211 FOR J=10 TO (NR-1@) STEP 190 

4212 PGXO=GX+I* PX*MX/ (NC=1) 

4213 PGX1=PGX6~-2*MX: PGX2=PGXO0+2*MX 

4214 PGYO@=GY+J* PY*MX/ (NR-1) 

4215 PGY1]1=PGYO@—-2*MX:PGY2=PGY0+2*MX 

4216 PRINT#2,"MA";PGX1;",";PGY@:PRINT#2, "DA"; PGX2;",";PGYO 

4217 PRINT#2,"MA";PGX8;",";PGY1:PRINT#2, "DA" ;PGXO;",";PGY2 


ee Ce Be Bo 


4218 
4228 
4236 
4245 
42590 
4252 
4254 
4256 
4258 
4268 
4261 
4262 
4263 
4264 
4270 
4271 
4272 
4273 
4274 
4275 
4288 


‘ 
4285 


4290 
4295 
4300 
4305 
4318 
4312 
4314 
4320 
4330 
4340 
4356 
4368 
4378 
4380 
4398 
4395 
4480 
4410 
4414 
4420 
4430 
4440 
4450 
4460 
4480 
4490 
4495 
4580 
4518 
4520 
TES 

4538 
454@ 
455¢@ 
4595 
4600 
4612@ 
4620 
46420 
46506 
4695 


NEXT J 
NEXT I 

GOTO 43886 

REM ***** Hewlett-Packard (HP-GL) ***#*# 
PRINT#2,°PS3;" 

PRINT#2,"IP 1780,602,15376,18602;" 
PRINT#2,"SC -10,188,-5,128;"% 
PRINT#2,°SP1;" 

ft] 


o 

;":TSEC=2 :GOSUB 60062 
":GOSUB 60806 

PRINT#2,"°PA";GX;",";PY;"3":GOSUB 68808 

PRINT#2,"PA";GX;","3GY;"3":GOSUB 68008 

FOR I=18 TO (NC-1@) STEP 18 

FOR J21@ TO (NR-10) STEP 18 

PGXO0=GX+I*PX/ (NC-1) 

PGX12PGXO0—2 s PGX2=PGX8+2 

PGYO=GY+J* PY/ (NR~1) 

PGY1=PGY@-2: PGY2=PGY6+2 


PRINT#2,"°PU;PA";PGX1;",";PGY6;";":PRINT#2,"PD; PA"; PGX2;","; PGYB 


PRINT#2,"PU;PA";PGXO0;",";PGY1;";":PRINT#2,"PD; PA"; PGX8;",";PGY2; 


NEXT J 

NEXT I 

ON HP GOTO 4700,4700,4318,4480 
REM ***#** Apple 410 **#es 

FOR J=1 TON 

IF SD%=1 THEN S(J) =P(J) 

IF SD%=2 THEN S(J) =Z (J) 

PRINT#2,"PS2;":PRINT#2,"LS2;" 


PRINT#2, "MA"; (GX+(X(J) —-XL) *MG) 3", "3 (GY+(Y¥(J)~YB) *MG) "3" 


PRINT#2,"PM8;”" 
PRINT#2,"LS1;" 
PRINT#2,"PL":PRINT#2,S (J) 
PRINT#2,CHRS$(13) CHRS$(18) CHRS$(3) 
NEXT J 
GOTO 4508 

REM ***** Epson HI-80 *****# 
FOR J=l1 TON 

IF SD%=1 THEN S(J)=P(J) 

YF SD%=2 THEN S(J) =Z(J) 
PRINT#2,"SP2":PRINT#2,"S140, 30" 


PRINT#2, "MA"; (GX+(X(J)~XL) *MG) *MX3","3 (GY+(Y(J) -YB) *MG) *MX 


PRINT#2,"AM8" 
PRINT#2,"SI25,20" 
PRINT#2,"LA";S(J) 
NEXT J 
GOTO 4608 

REM ***#* Anple 410 **#e® 
PRINT#2,"LS2;" 

PRINT#2,"MA 10,-~-3;" 


PRINT#2,"PL":PRINT#2,"PC-Random/RData File ~— ";FS;" ; Date 


PRINT #2,CHRS$(13) CHRS$(18) CHRS(3) 
CLOSE#2 
GOTO 200¢ 

REM ***** Epson HI-80 ***#*® 
PRINT#2, "S146, 30" 
PRINT#2,"MA";5*®MX;","3 2*MX 


PRINT#2,"LAPC-Random/RData File — ";FS$;" : Date :";DATES 


CLOSE#2 
GOTO 2000 
REM ***** Hewlett-Packard (HP-GL) ***#*# 


4760 
4712 
4714 
4726 
4730 
4750 
4760 
4770 
4780 
4800 
4818 
4820 
TES 

4830 
4840 
4858 
4995 
5995 
6000 
6016 
6028 
6030 
6048 
6168 
6128 


FOR J=1 TON 

IF SD%=] THEN S(J) #=P(J) 

IF SD@=2 THEN S(J)#Z(J) 

PRINT#2,"SP3;" 

PRINT#2,"PU; PA"; (GX+(X(J) -XL) *MG) 5","% (GY+(¥(J)~YB) *MG) 7"; 
PRINT@2,°SI0.18,0.22;" 

PRINT$2,"LB+";S(J) 

PRINT#2,CHRS (3) 

NEXT J 

PRINT#2,"SP2;S10.2,0.3;" 

PRINT#2,"PU;PA 18,-3;" 
PRINT#2,"LB":PRINT#2,"PC~Random/RData File -— “;F$;" : Date :";DA 


PRINT#2,CHRS (3) 
CLOSE#2 
GOTO 2008 
REM ***8* User DATA Storage area ***#* 
REM ***## Printout of Data File Details ***** 
CLS 
GOSUB 9860 
PRINT ° Place paper in Printer and switch on” 
PRINT 
INPUT " To Continue Press <Enter> ";R$ 


LPRINT "PC-Random/RData File — ";F$;" - Date :";DATES 
LPRINT QRRARARRAARARAAKARARAARAARAARARARAAARARARARARARARARARARKARRARAERKRAERKAARAKAALR 


thanhanhanannee 


6130 
6146 
6150 
6166 
6176 
6188 
6198 
6208 


LPRINT 

LPRINT "X=-coordinate Left hand side ";TAB(35) ;" 

LPRINT "“X=coordinate Right hand side";TAB(35) ;”* 
as 
4" 


LPRINT “Y-coordinate at foot of plot";TAB(35 = 
LPRINT °Y-coordinate at top of plot";TAB(35 = 
LPRINT “Number of data points”";TAB(35) ;"=";N 
LPRINT 


LPRINT "Point No.";TAB(15) ;"X-coord”;TAB(36) ;"Y-coord”;TAB(45) ;"* 


Z2-coord"” 


6218 
6226 
6230 
6246 
6258 
6268 
6270 
8995 
98600 
9818 
9028 
9838 
9046 
9858 
9588 
9582 
9564 
9506 
9508 
9518 
9520 
9530 
9548 
9550 
9560 
9570 
9588 
9598 
9608 
9658 


CPRINT “2282s sSS2 SSS StS SSS ees eet eee SSS St es Sess sssssszsssssseee==" 
FOR I=] TON | 

LPRINT P(I) ;TAB(14);X(I) ;TAB(29) ;Y(1) ;TAB(44) 7Z(1) 
nee SHAR AAARAARAAAAAARARAAAARHARRAARAKRARARRARKRKAARARKAKRAKKARKRAKKKKRaAKSE! 
LPRINT 

SR%=1:GOTO 2608 

REM ***** pC-Random Heading ***** 

PS="PC-Random : ":TS=PS$+F$ 

LOCATE 1,2:PRINT TS 

LOCATE 1,30:PRINT DATES 

PRINT ® seseeee2e22222222282222222222e2e2ee2e2s===" 

PRINT 

RETURN 

CLS 

SCREEN SCR%:IF SCR%=@0 THEN 9518 

IF SCR%=1]1 THEN 9568 

SCREEN 1:GOTO 9518 

COLOR 8,90 

GOSUB 98890 

PRINT:PRINT * Program Complete ” 

PRINT 

PRINT:PRINT ™ Do you wish to run PC-Contour now (y/n) "® 
INPUT QS:IF QS$="n" THEN 9680 

IF Q$<>"y" THEN 9540 

PRINT: PRINT:PRINT “ Loading PC-Contour from disc .....” 
RUN "CONTOUR® 

PRINT 

PRINT: PRINT 

END 
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3.3.5 PC-Random/Demo 


To be able to run PC-Random/RData it is necessary to prepare a data 
file. It is suggested that once the following DATA statements have been 
entered into the file PC-Random/RData that it is resaved under the title 
“RANDDEMO” with changed remarks in lines 10 and 25. The data file 
presented is for a tacheometric survey at Glenmore, and the resulting data 
file will be saved automatically on disc as “GLENMORE” for subsequent 
access by PC-Contour and PC-Surfplot. 

Initial data statements give the file name “Glenmore” followed by the 
X-min, X-max, Y-min, Y-max coordinates of the plotted area. This is then 
followed by the number of data points (30), and the X, Y, Z-coordinates of 
each point. For simplicity, it is suggested that each data point 1s entered ina 
separate DATA statement for ease of checking and to save confusion with 
the data points. There is no need to identify the maximum and minimum 
Z-values as these will be located in the READ routine. 


3.3.6 PC-Random/Demo -— data 


5000 DATA "Glenmore" 

5818 DATA 6,400 

5028 DATA 8,380 

5030 DATA 58 

584@ DATA 1,372,24,173.3:REM Line 1 
505@ DATA 2,3708,74,176.8 

5068 DATA 3,368,122,185.6 

5078 DATA 4,366,148,187.7 

5088 DATA 5,364,185,192.9 

5098 DATA 6,362,218,204.2 

5188 DATA 7,360,254,213.8 

5118 DATA 8,358,273,216.8 

5128 DATA 9,342,97,177.2:REM Line 2 
5138 DATA 16,309,122,178.5 

51408 DATA 11,272,148,188.4 

5158 DATA 12,241,171,182.6 

5168 DATA 13,205,196,186.5 

5178 DATA 14,176,219,189.3 

5188 DATA 15,144,242,194.8 

5198 DATA 16,121,261,204.1 

5208 DATA 17,99,278,208.3 

5210 DATA 18,251,276,209.2 

5220 DATA 19,336,41,174.3:REM Line 3 
5238 DATA 20,304,57,188.9 

524@ DATA 21,269,76,186.1 

$258 DATA 22,228,101,188.8 

5260 DATA 23,205,122,191.3 

5278 DATA 24,172,131,193.7 

5288 DATA 25,147,143,196.2 

5298 DATA 26,169,166,199.1 

5388 DATA 27,72,185,202.4 

5318 DATA 28,61,216,265.8 

5328 DATA 29,45,260,212.2 

5338 DATA 30,205,32,178.5 

53348 DATA 31,236,41,182.3:REM Line 4 
5358 DATA 32,287,101,183.6 
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5368 DATA 33,339,156,186.5 

5378 DATA 34,133,13,178.3:REM Line 5 
5388 DATA 35,155,42,182.8 

5398 DATA 36,182,86,188.1 

5400 DATA 37,224,136,187.2 

5418 DATA 38,265,196,189.3 

5420 DATA 39,288,229,2068.9 

5438 DATA 48,369,259,211.1 

5448 DATA 41,329,277,215.7 

5458 DATA 42,21,26,178.2:REM Line 6 
54608 DATA 43,54,88,187.3 

5478 DATA 44,76,119,193.9 

5488 DATA 45,132, 282,193.8 

5498 DATA 46,171,256,197.4 

5588 DATA 47,187,277,286.2 

5518 DATA 48,96,51,182.9 

5528 DATA 49,135,116,193.5 

5538 DATA 568,229,238,196.2 


3.3.7 PC-Random — computer printout of Demo file 


PC-Random/RData File -—- Glenmore - Date :981-81-1987 
RARAAAARAARARAAARRARAARAAAARARARAAARAKARRARARARARAARARARARARAKAKEAKEERARARARKEKAK 


X-coordinate Left hand side = @Q 

X-coordinate Right hand side = 480 

Y-coordinate at foot of plot = 6 

Y-coordinate at top of plot = 300 

Number of data points = 56 

Point No. X-coord Y-coord zZ-coord 
I 372 24 173.3 

370 74 176.8 

3 368 122 185.6 
4 366 148 187.7 
5 364 185 192.9 
6 362 218 204.2 
7 360 254 213.8 
8 358 273 216.8 
9 342 97 177.2 
1@ 309 122 17825 
ll 272 148 188.4 
12 241 171 182.6 
13 265 196 186.5 
14 176 219 189.3 
is 144 242 194.8 
16 121 261 204.1 
17 99 278 288.3 
18 251 276 269.2 
19 336 41 174.3 
208 304 57 1898.9 
21 269 76 186.1 
22 228 161 188.8 
23 205 122 191.3 
24 172 131 193.7 
25 147 143 196.2 
26 189 166 199.1 
27 72 185 282.4 
28 61 216 205.8 
29 45 266 212.2 
36 2065 32 178.5 


31 236 41 182.3 
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32° 287 101 183.6 
3a 339 156 186.5 
34 133 13 178.3 
35 155 42 182.8 
36 182 88 — 188.1 
37 224 136 187.2 
38 265 196 189.3 
39 288 . 229 200.9 
40 389 259 211.1 
41 329 277 215.7 
42 21 26 178.2 
43 54 88 187.3 
44 76 119 193.9 
45 132 202 193.8 
46 171 256 197.4 
47 187 277 206.2 
48 96 51 182.9 
49 135 118 193.5 
58 229 238 196.2 


RARARAARARARAAARAKRARARRAKARKRERRARARKRARKEKARARARARAKKEKKEKRARKREAR 


3.3.8 PC-Random/FData 


This is the third program in the PC-Random suite and is provided for those 
who may already have random data stored in a disc file, for example, as 
generated by two programs in the PC-Survey suite of programs for land 
surveying. Two of the programs from BASIC Programs for Land Survey- 
ing allow field data collected by land surveying or topographic mapping, 
using either tacheometry or EDM (Electromagnetic Distance Measure- 
ment) to be plotted on various plotters and stored on disc for access by 
PC-Random/FData. 

This program follows the same format as PC-Random/R Data and only the 
lines listed in Section 3.3.9 need to be altered. The main changes occur in 
lines 1010-1060, where instead of reading data statements, the data is read 
from a disc file. 


3.3.9 PC-Random/FData — BASIC program 
Changes to PC-Random/RData 


18 REM <RANDFDAT> Program for IBM PC and Compatibles 
25 REM Reads DATA File stored under FS+".DAT" by PC-Survey 


178 PRINT " Select application required :-" 


188 PRINT 

198 PRINT " <1> Read DATA and present details® 
288 PRINT 

218 PRINT " <2> Read DATA and store on disc® 
228 PRINT 


238 PRINT " <3> Quit °; 
268 INPUT RAND@: IF RAND&<1 OR RAND%>3 THEN 160 
278 IF RAND&=3 THEN 9580 
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995 REM **### DATA File Read Routine **#*#*# 


1880 
1018 
1028 
1030 
1048 
1058 
1068 


1130 
11480 
1158 
1168 
1165 
3218 


4520 
TES 


46208 


48298 
TES 


6180 


CLS:GOSUB 9888 

INPUT ® Enter DATA file name °;FS 
D$S#FS$+".DAT® 

OPEN °"1I°,#4,D$ 

INPUT#4,XL,XR 

INPUT#4,YB,YT 

INPUT#4,N 


INPUT#4,P(I),X(I),Y¥(1I),2(1) 


IF Z(I) > MAX THEN MAX = 2(1I) 
IF Z(I) < MIN THEN MIN = Z(TI) 
NEXT I 
CLOSE#4 
LOCATE 1,2:PRINT "PC-Random/FData — °;FS$ 


PRINT#2,"PL":PRINT#2,"PC-Random/FData File — ";F$;" : Date : 


PRINT#2,"LAPC-Random/FData File -— ";F$;" : Date :";DATES$ 


PRINT#2,"LB”":PRINT#2,"PC-Random/FData File — ";F$;" : Date : 


LPRINT “PC-Random/FData File -— ";FS$;" —~ Date :";DATES$ 


3.4 PC-CONTOUR 


This program accesses disc files generated by either PC-Grid or PC-Random 
and can be run with the two demonstration data files supplied, that 1s 
“GLAPWELL” for spoil-heap contours and “GLENMORE” for a 
tacheometnic survey. 

On running PC-Contour, the user is asked to enter the file name for the 
disc data file. After reading the data file, the data information is displayed on 
screen with advice on the size of the plot, the numbers of rows and columns 
in the grid, and the maximum and minimum values for the plot. These latter 
values will be the limits for the contour values when plotting on screen or 
plotter. 


3.4.1 PC-Contour — subroutine index 


(a) 
(b) 
(c) 
(d) 
(e) 
(f) 
(g) 


Line 
numbers Function 
10—90 Initialization and control 


100-1170 Data read routine 
2000-2180 Data display and plot selection 
2500-3990 Graphics screen display 
4000-5190 Plotter routines 
9000-9050 Routine to clear screen 
9500-9540 Termination of program 
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(a) Initialization and control Line numbers 10-90 


The program is initialized and the installation file checked for the availability 
of graphics screens and colour for subsequent display of the contours. 


(b) Data read routine - Line numbers 100-1170 


This routine requests the data file name and prompts the user with the two 
demonstration data files prepared in Sections 3.2 and 3.3. 


(c) Data display and plot selection Line numbers 2000-2180 


The presentation of the data information is similar to that described in 
Sections 3.2 and 3.3. The user should note the maximum and minimum 
Z-values to be used in the selection of contour levels, before selecting either 
the graphics screen or plotter for the contour display. 


(d) Graphics screen display Line numbers 2500-3990 


The user is first given a choice of low-resolution (SCREEN 1) or 
high-resolution (SCREEN 2) for the graphics display which uses the format 
shown in Fig. 1.4. In addition to the normal box outline, grid nodes are 
plotted at the junction of every ten grid lines. 

The right-hand text column ts used to enter the lowest contour level, the 
highest contour level and the contour interval. Contours are interpolated 
between each set of four grid points in turn using an inverse square distance 
weighted interpolation. During plotting, the current contour value is 
displayed on screen. Each time a new contour is encountered its value will be 
noted alongside the contour. It is possible at the end of the plot to go back 
and insert more contours or to start again with different contour values. 

Thus it is suggested that for the first plot, a contour interval of 10 is 
selected, followed by 5, for example. As the graphics plotter routines take 
longer than the screen plot, it is often useful to run the screen display first to 
check the output before using the plotter. On completion, a screen dump of 
the graphics display can be obtained on a dot-matrix graphics printer 
(Fig. 3.3) for comparison with the data (Figs 3.1 and 3.2) of the PC-Random 
file “GLENMORE”. The PC-Grid plot file “GLAPWELL” is shown on 
Fig. 3.4. 


(e) Plotter routines Line numbers 4000-5190 


Three plotter routines are provided for the Apple 410, Epson HI-80 and 
Hewlett-Packard HP-7475A colour plotters as discussed in Chapter 2. Most 
of the colour plotters take four pens and the default settings are 1-black, 
2-red, 3-green, 4-blue. The border frame and grid nodes are plotted by pen 1 
(black), the 10-unit contour is plotted by pen 3 (green), the 5-unit contour by 
pen 4 (blue) and the 2.5-unit contour by pen 2 (red). Each time a new 
contour is encountered its values will be noted alongside the contour. All 
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Fig. 3.3. PC-Contour plot of contours for Glenmore tacheometry survey on Screen 1 
(medium-resolution). 


other unit values (e.g. 2) will be plotted by pen 3 (green). The final title, 
location and date will be plotted by pen 2 (red). If, however, the user wishes 
to change these colours, the appropriate pen selection should be changed 
using the commands in Table 2.1. 

The contour plot for the PC-Random file “GLENMORE” is shown in 
Fig. 3.5 for comparison with Fig. 3.3. The contour plot for the PC-Grid 
file “GLAPWELL” is shown in Fig. 3.6 for comparison with Fig. 3.4. 


(f) Routine to clear screen Line numbers 9000-9050 


This routine is provided to clear the screen before each plot and provide text 
at the top of the screen to identify the program title and data file presentation 
together with the current date. 


(g) Termination of program Line numbers 9500-9540 


This is the final display at the completion of the program and the user is given 
an option to run PC-Contour again or return to BASIC command level. 
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Fig. 3.4 PC-Contour plot of Glapwell spoil-heap contours of Screen 2 (high- 
resolution). 


3.4.2 PC-Contour — numeric variables 


A = Variable used in contour interpolation 

B Vaniable used in contour interpolation 

C Variable used in contour interpolation 

CC Current contour value 

CCOL Column number for text display on graphics screen 
CROW Row number for text display on graphics screen 


Variables used in contour interpolation 


~~) 
r -) 
Ie to 


DX X-distance of plotted area 

DY = Y-distance of plotted area 

GPX98-2 = Vaniables to plot grid nodes on screen 
GPY8-2 = Vanables to plot grid nodes on screen 

Gx = X-coordinate at lower left-hand corner of plot 
GY = Y-coordinate at lower left-hand corner of plot 
HC = Highest contour on plot 

HP% = Integer variable to identify plotter 

H(NR,NC) = Array to read Z-values 
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Fig. 3.5 PC-Contour plot of contours for Glenmore tacheometry survey on plotter. 


A 


LC 
MAX 
MG 

MH 
MIN 
MV 

MX 

NC 

NP 

NR 
PGXH—2 
PGY#-—2 
PP 

PX 

PY 
PLOT % 
RGB&% 
SN 

SY 
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Integer variables used in loops 


Lowest contour on plot 


= Maximum Z-value in data 


Multiplication factor for plotted display 
Multiplication factor for honzontal resolution 
Minimum Z-value in data 

Multiplication factor for vertical resolution 
Multiplication factor for Epson HI-80 plotter 
Number of columns tn grid 

Number of points in contour interpolation 
Number of rows in grid 

Variables to plot grid nodes on plotter 
Variables to plot grid nodes on plotter 
Screen display variable 

X-coordinate at upper right-hand corner of plot 


Y-coordinate at upper right-hand corner of plot . 


Integer variable for plotter selection 
Integer variable to select line colour 
Horizontal display width in user units 


= Vertical displav height in user units 


Fig. 3.6 


SCD% 
SCR& 
SR%& - 
TL 
X,X9,X1 
Y,Y6,Y1 
Z1-Z4 
ZT 
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PC-Contour plot of Glapwell spoil-heap contours on plotter. 


= Integer variable to define screen display 

= Integer variable to select screen or plotter 

= Integer variable for screen resolution 

= Variable used in titling contour level 

= Variables used in contour interpolation 

= Variables used in contour interpolation 

= Z-levels of grid points in contour interpolation 
= Variable used in contour interpolation 


3.4.3. PC-Contour — string variables 





DATES 


F$ 
GD$ 
M$ 
P$ 
PCF$ 
Q$ 
T$ 


Variable used to return the date of the computer's calendar/clock in 
US format MM-DD-YYYY. 
File name 
Entry for graphics dump on pmnter 
= “MORE”, “ERASE” or “END” 
= “PC-Contour”’ 
= “PCSCREEN”’, installation file for graphics 
Question response (Y/N) or (y/n) 
= PS + FS 





LUO SURFACE MODELLING 


3.4.4 PC-Contour — BASIC program 


10 REM <CONTOUR> Program for IBM PC and Compatibles 
12 REM Version PC-1.8, @1-01-1987 

14 REM (C) Copyright 1987 P.H. Milne 

16 REM All rights reserved 

2@ REM CONTOUR Plot on <1l> Apple 410, <2> Epson HI-88, or <3> HPGL 
Colour Plotter 

22 REM Setting in line 4898 "%COM1:1200 N,7,1,CS2008,DS15000" 
25 REM This version includes Graphics printer dump of Screen 
3@ CLEAR 

35 PCFS$="PCSCREEN® 

4@ OPEN "I°,#3,PCFS$ 

45 INPUT#3,SCR& 

5@ INPUT#@3,SCD% 

55 CLOSE#3 

6@ ON (SCR%¥+1) GOTO 62,64,66 

62 SCREEN @:WIDTH 408:GOTO 80 

64 SCREEN 1:GOTO 68 

66 SCREEN 1:GOTO 88 

68 IF SCD%&=]1 THEN 828 

78 COLOR 8,8 

8@ KEY OFF 

90 FS=" ";RGB%=0:SR¢=) 

18@ CLS 

11@ GOSUB 9800 

128 PRINT 

138 PRINT *® THIS PROGRAM READS AN ELEVATION FILE® 
148 PRINT " FROM A REGULAR SQUARE GRID TO DISPLAY” 
15@ PRINT “ A CONTOUR MAP ON SCREEN OR PLOTTER." 

168 PRINT 

268 GOSUB 1000: REM Read DATA FILE 

218 GOTO 28008: REM Select Output 

995 REM ***** Read Data from FILE ***** 

1808 PRINT:PRINT ® PC-Grid/Demo file called GLAPWELL":PRINT “" PC-R 
andom/Demo file called GLENMORE” 

1005 PRINT: PRINT:PRINT “ Enter DATA FILE name ‘"; 
1810 INPUT FS:IF FS="" THEN 10805 

1015 PRINT: PRINT:PZINT " One moment please ..... reading data":PRINT 
1020 OPEN "I",#1,FS 

1038 INPUT#1, NR 

18048 INPUT#1l, NC 

1058 INPUT#1l, MAK 

1860 INPUT#1, MI?. 

187@ INPUT#1, DX 

1088 INPUT#1, CY 

1098 DIM H((NR - i),(NC - 1)) 

110@ FOR I = (NR-:) TO @ STEP -1 

1110 FOR J = @ TO (NC - 1) 

1120 INPUT#l, Hl 

1130 H(I,J) = Hl 

11460 NEXT J 

1158 NEXT I 

1160 CLOSE #1 

1170 RETURN 

1995 REM ***** Data Information and Plot Selection ***** 
2006 CLS 

2802 SCREEN SCR%:IF SCR*=8 THEN 2012 

2004 IF SCRR=1]1 THES 2068 

2006 SCREEN 1:GOTC 2018 

2088 COLOR 8,8 

2018 GOSUB 9088 

2028 PRINT “ THE DATA FILE °;F$;" HAS BEEN READ" 
280308 PRINT 


2848 PRINT ® No. Rows = ";NR;" No. Cols.= “;NC 
2058 PRINT ” Minimum value = ";MIN 
286@ PRINT " Maximum value = ";MAX 


2870 
20880 
28690 
21900 
2118 
2128 
2130 
2140 
2158 
2168 
2172 
2188 
2495 
2508 
2510 
2528 
2530 
2548 
2558 
2568 
2570 
2588 
2608 
2995 
3008 
3010 
38208 
3630 
3648 
3858 
3860 
3870 
30828 
3098 
3100 
3112 
3128 
3130 
3140 
3158 
31690 
3170 
3182 
3198 
3195 
32080 
3220 
3230 
3240 
3250 
3252 
3255 
3260 
3262 
3265 
3270 
3272 
3275 
3288 
3290 
3320 
3310 
3312 
3315 
3328 
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PRINT " Range = "sINT( (MAX — MIN) *100+.1)/108 
PRINT " 1 Map row = ";(DY / (NR-1));" UNITS® 
PRINT * 1 Map col. = "3;(DX / (NC-1));" UNITS” 

BEEP 


PRINT: PRINT:PRINT " Select Display :-" 

PRINT:PRINT " <1> SCREEN, <2> PLOTTER or <3> QUIT” 
IF SCR&<> @ THEN 2158 

PRINT: PRINT ® N.B. Option <l> not available":PRINT 
INPUT SELECTION 

IF SELECTION<1 OR SELECTION>3 THEN 2108 

IF SCR*=0 THEN IF SELECTION=1 THEN 2000 

ON SELECTION GOTO 2562,4808,9500 


REM Choose Screen for Plotting 

CLS 

GOSUB 9880 

PRINT:PRINT TAB(3) ;"Select SCREEN";TAB(28);"<l> Colour —-—" 
PRINT TAB(2@);" Low resolution" 

PRINT: PRINT TAB(28);"<2> Black & White -" 

PRINT TAB(2@) ;” High resolution” 


PRINT TAB(3);: INPUT SR&% 
IF SR&¥<] OR SR¥>2 THEN 2529 


IF SR*¥=1 THEN 3000 

SCREEN 2 

REM ***** Screen Plot ***** 
CLS 


GX=5*SRB:GY=5 sMG=0 : MH=6 :MV=0 :SX=240*SR%:SY=1 80: PP=0 
MH=SX/ ((NC~1) *SR&) :MV=SY/ (NR=-1) 

IF MH<=MV THEN MG=MH ELSE MG=MV 
WINDOW (0,8) = (320*SR%, 200) 

IF SR%=]1 THEN RGB%=2 ELSE RGB&%=3 
PX=GX+ (MG* (NC~1) *SR&) : PY=GY+MG* (NR=]) 
LINE (GX,GY)-(PX,PY) ,RGB&,B 

IF SR%=1 THEN 3102 

LINE (GX-1,GY)~-(PX-1,PY) ,RGB&,B 

FOR I=18 TO (NC-1@8) STEP 18 

FOR J=10 TO (NR-16) STEP 1@ 

GPX8=GX+ (MG*I*SR8) 

GPX1=GPXO-5*SR& :GPX2=GPX0+5*SRE 


GPY@=GY+(MG*J) 
GPY1=GPY@-S5:GPY2=GPY0+5 

LINE (GPX1,GPYG6) - (GPX2,GPY@) ,RGB% 
LINE (GPX@,GPY1) ~ (GPX2,GPY2) ,RGB% 
NEXT J 

NEXT I 


IF PP=1. THEN 3888 

" Start of ‘MORE’ contours 
LOCATE 1,2:PRINT "PC-Contour = ";FS 
LOCATE 1,30*SR%&:PRINT DATES 

LOCATE 3,32*SR&:PRINT “CONTOURS” 
LOCATE 4,32*SR&:PRINT "LOW VAL:" 
LOCATE 5,32*SR&:PRINT ” - 
LOCATE 5,32*SR&: INPUT LC 

LOCATE 6,32*SR&:PRINT °HIGH VAL:” 
LOCATE 7,32*SR&:PRINT * " 
LOCATE 7,32*SR%&: INPUT KC 

LOCATE 9,32*SR%:PRINT "CONTOUR" 
LOCATE 180,32*SR%&: PRINT *% = 
LOCATE 10,32*SR*: INPUT “INT.";CI 
IF LC < (MIN - CI) THEN LC = LC + CI 
IF HC>MAX THEN HC=MAX 

FOR CC = LC TO HC STEP CI 

LOCATE 11,32*SR@: PRINT “CURRENT” 
LOCATE 12,32*SR&:PRINT ” " 
LOCATE 12,32*SR&:PRINT "=";CC 

FOR I = @ TO (NR-2) 
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3330 
3348 
3358 
33690 
3372 
3380 
3398 
3400 
3418 
3420 
3430 
3440 
3450 
3460 
3470 
3488 
3490 
3586 
3510 
3522 
3538 
3548 
3558 
3568 
3570 
3580 
3598 
3600 
3618 
3620 
3638 
3631 
3632 
3633 
3634 
3635 
3636 
3637 
3638 
3639 
3646 
3650 
3660 
3670 
3680 
3690 
3708 
3710 
3720 
3736 
3748 
3758 
3760 
3778 
3788 
3798 
3800 
3818 
3820 
3830 
3848 
3858 
3868 
387¢ 
3875 
3880 
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Y@ = MG* I 
FOR J = @ TO (NC - 2) 
X@ = MG* J 
NP = @ 
Z1 = H(I,J) 
IF CC > Z1 THEN NP = NP + 1 
Z2 = H(I, (J + 1)) 
IF CC > Z2 THEN NP = NP + 1 
Z3 = H((I + 1),J) 
IF CC > Z3 THEN NP = NP + 1] 
24 = H((I + 1),(J + 1)) 
IF CC> 24 THEN NP = NP + 1] 
IF NP = @ OR NP = 4 GOTO 3848 
A = 2] 
B= 22-A 
C = 23 -A 
D=Z4-+A-B-+C 
ZT = O 
FOR Yl = @ TO 1 STEP .25 
DR = B + D* Yi 
IF DR = @ GOTO 36828 
Xl = (CC - A =-C * Y¥1) / DR 
IF Xl < @ OR Xl > 1 GOTO 36828 
X = X86 + MG* Xl 
Y = YO + MG* Yl 
IF X > PX OR Y > PY GOTO 36790 
IF ZT = @ THEN 3620 
LINE - ((X*SR%+GX) ,(Y+GY)),1 
GOTO 36786 


PSET ((X*SR%+GX) , (Y+GY)),1 

IF TL>@ THEN 3670 

IF CC/18 <> INT(CC/16) THEN 3670 
CCOL= (X*SR&+GX) /8: CROW=(Y+GY) /8 

IF CCOL<2*SR& THEN CCOL=2*SR% 

IF CCOL>27*SR%& THEN CCOL=27*SR8% 

IF CROW<=2 THEN CROW=2 

IF CROW>22 THEN CROW=22 

IF CC/28 <> INT(CC/28) THEN 3639 
LOCATE (25-CROW) ,CCOL:GOTO 3642 


LOCATE (24-CROW) ,CCOL 
PRINT CC 
TL=TL+1 
GOTO 3620 
ZT = 2T + 1 
NEXT Yl 
ZT = @ 
FOR Xl = @ TO 1 STEP .25 
DS = C + D *® Xl 
IF DS = 86 GOTO 3839 
Yl = (CC - A - B * X1) / DS 


IF Yl < @ OR Y1 > 1 GOTO 3838 
X = X08 + MG* X1 
Y = Y¥@. + MG* Yi 
IF X > PX OR Y > PY GOTO 38290 
IF ZT = @ THEN 3819 
LINE - ((X*SR&+GX), (Y+GY)),1 


GOTO 38220 
PSET ((X*SR8+GX) ,(Y+GY)),1 
ZT = ZT + 1 
NEXT Xl 
NEXT J 
NEXT I 
TL=0 
NEXT CC 
PP=1:GOTO 30676 


LOCATE 14,32*SR&@:PRINT "SCREEN" 


3882 
3884 
3885 
3886 
3888 
3898 
3960 
3910 
3928 
3938 
3948 
3945 
3950 
3955 
3960 
3970 
3986 
3990 
3995 
4000 
4810 
4828 
4830 
4040 
4058 
4051 
4852 


4054 
4056 
4658 
40859 
4670 
40880 
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LOCATE 15,32*SR¢:PRINT "“DUMP(y/n) ” 
LOCATE 16,32*SR&:PRINT * " 
BEEP 
LOCATE 16,32*SR&: INPUT GDS$ 
IF GD$="n" THEN 3918 
IF GDS$<>"y" THEN 3888 
LCOPY 
LOCATE 18,32*SR&:PRINT “™"MORE'" 
LOCATE 19,32*SR%&:PRINT “CONTOURS” 
LOCATE 28,32*SR&:PRINT “*ERASE‘," 
LOCATE 21,32*SR&:PRINT “or ‘END‘" 
BEEP 
LOCATE 22,32*SR%&:PRINT " ” 
LOCATE 22,32*SR®&: INPUT MS 
IF MS = “MORE" GOTO 328690 
IF MS = "ERASE" GOTO 3008 
MS=" ";SR%=1 
CLS:GOTO 2000 
REM ***** Colour Plotter Display ***** 
IF MS="MORE" THEN 4238 
GX=@0:GY=0 :MG=0:MH=@ : MV=6 :SX=192:SY=1 20 
MH=SX/ (NC-1) :MV=SY/ (NR-1) :MX=2670/SX 
IF MH<=MV THEN MG=MH ELSE MG=MV 
PX=MG*® (NC-]1]) : PY=MG*® (NR-1) 
CLS: GOSUB 9800 
PRINT: PRINT ”® Select Plotter <l> Apple 416" 
PRINT:PRINT ” <2> Epson HI-80":PRINT:PRINT " 
<3> Hewlett-Packard “;:INPUT PLOTS 
IF PLOT%<l OR PLOT%>3 THEN 4050 
IF PLOT%=1 THEN HP=3:REM Apple 
IF PLOT%=2 THEN HP=4:REM Epson 
IF PLOT%=3 THEN HP=1:REM Hewlett-Packard (HPGL) 
PRINT: PRINT:PRINT " Place paper on Plotter and switch on" 
PRINT: PRINT " To Continue Press <Enter> ":WHILE RS="":RS=INKEYS 


:WEND 


4896 
4100 
4163 
41¢@5 
4118 
4122 
4148 
4141 
4142 
4143 
4144 
4150 
4155 
4160 
4161 
4162 
4163 
4178 
4175 
4180 
4185 
4190 
4195 
4200 
4201 
42082 
4203 
4204 
4205 


OPEN "COM1:120@ N,7,1,CS2000,DS15000" FOR OUTPUT AS #2 

ON HP GOTO 5108,51080,4105,4200 

REM ***** Apple 4)Q0 **#t# 

PRINT#2,"SP3;” 

PRINT#2,"WD -5,-5,192,125;" 

PRINT#2,"PS1;" 

PRINT#2,"MA":GxX;",";GY;3";" 
y FON" e* 

PRINT#2,"DA";:PX;",";PY;";" 
f° SPYs" =™ 


FOR I=18 TO (NC-18) STEP 198 

FOR J=10@ TO (NR-18) STEP 18 

PGXO=GX+ (MG*#I) 

PGX1=PGXO6-2: PGX2=PGX8+2 

PGY@=GY+(MG*J) 

PGY1=PGY@-2: PGY2=PGY8+2 
PRINT#2,"MA";PGX1;",";PGYG;";":PRINT#2, "DA"; PGX2;","%;PGY0;";" 
PRINT#2,"MA";PGXO;",";PGY1;";":PRINT#2,"DA";PGXO;",";PGY2;";" 
NEXT J 

NEXT I 

GOTO 42390 

REM ***** Epson HI-8@ ****# 

PRINT#2,"IWG,8,2670,1926" 

PRINT#2, "MA" ;GX;",";GY 

PRINT#2, "OR" 

PRINT#2,"SP1":PRINT#2, "LTO" 

PRINT#2,"MA";GX;",";GY 
PRINT#2,"DA";PX*®*MX;"7,"sGY3",":PX*MX;","sPY*MX:","%:GX3","3;PY*MX;"% 


0" 3GX3","sGY 


4210 


FOR I=16 TO (NC-16) STEP 18 


SURFACE MODELLING 


FOR J=18 TO (NR-1@0) STEP 18 
PGX@=GX+ (MG*I *MX) 
PGX1=PGXO0-2*MX:PGX2=PGX0+2*MX 
PGYO=GY+(MG*J*MX) 
PGY]=PGY@-2*MX: PGY2=PGY0+2*MX 
PRINT@2, "MA"; PGX1;",";PGYO:PRINT#2, "DA"; PGX2;","; PGYB 
PRINT#2,"MA";PGX0;",";PGY1:PRINT#2, "DA"; PGXO0;",";PGY2 
NEXT J 
NEXT I 
REM ==2s2= Return to SCREEN ===== 
CLS:GOSUB 9800 
LOCATE 8,10:PRINT “CONTOUR LIMITS:”* 
LOCATE 10,10:INPUT " LOW VALUE "“;LC 
LOCATE 12,10: INPUT " HIGH VALUE ";HC 
LOCATE 14,1@:INPUT "CONTOUR INTERVAL =";CI 
IF LC < (MIN - CI) THEN LC = LC #+ CI 
IF HC>MAX THEN HC=MAX 
FOR CC # LC TO HC STEP CI 
LOCATE 16,10:PRINT “CURRENT CONTOUR =";CC 
IF CC/10=INT(CC/10) THEN 4305 
IF (CC-5)/10=INT((CC-5)/18) THEN 4310 
IF (CC-2.5) /S=INT( (CC-2.5)/5) THEN 4315 
ON HP GOTO 4306,4306,4308,4309 
PRINT#2,"SP3;":PRINT#2,"LT;":GOTO 43220 
PRINT#@2,"PS3;":PRINT#2,"LT1;"%:GOTO 4320 
PRINT#2,"SP3":PRINT#2,"LTO":GOTO 4320 
ON HP GOTO 4311,4311,4313,4314 
PRINT#@2,"SP4;":PRINT#2,"LT;":GOTO 4320 
PRINT#@2,"PS4;":PRINT#2,"°LT1;":GOTO 4328 
PRINT#2,"SP4":PRINT#2,"LTO":GOTO 4320 
ON HP GOTO 4316,4316,4318,4319 
PRINT#@2,"SP2;":PRINT#2,"LT;":GOTO 43280 
PRINT#2,"PS2;":PRINT&#2,"LT1;":GOTO 4320 
PRINT#2,"SP2":PRINT#2,"“LTO" 
FOR I = @ TO (NR-2) 


YO = MG* I 
FOR J = @ TO (NC - 2) 
XO = MG* J 
NP = oO 
Zl 2 H(I,J) 


IF CC > 21 THEN NP 
£2 = Hitetd + 1)) 

IF CC > 22 THEN NP = NP + 1 
23 = H((T + 1),2) 

IF CC > 23 THEN NP = NP + 1 
z4 = H((I + 1),(3 + 1)) 

IF CC> 24 THEN NP = NP + 1 

IF NP = @ OR NP = 4 GOTO 4848 


NP + 1] 


A = 21 

B=2Z2-A 

C = 23-=] A 
D=Z24-A-8B-C 
ZT = 2 


FOR Yl = @ TO 1 STEP .25 
DR =B+D* Yi 

IF CR = B@ GOTO 46890 
Xl = (CC ~ A = C * Y1) / DR 

IF Xi < 6 OR Xl > 1 GOTO 4688 
X = Xd + MG* Xl 
Y = Ye + MG* Yi 

IF X > PX OR Y > PY GOTO 46728 

IF ZT = @ THEN 4615 
ON HP GOTO 46@01,46H01,4603,4605 
PRINT#2,"PD;PA";X+GX3","3Y+GY;";":sGOTO 4670 
PRINT#@2,"DA";X+GX;","3Y+GY3;";":GOTO 4670 
PRINT#2,"DA"; (X+GX) *MX3","s (Y+GY) *MX 
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4618 GOTO 4678 

4615 ON HP GOTO 4616,4616,4618,4619 

4616 PRINT#2,"PU;PA";X+GX;",";Y+GY;";":GOTO 4620 
4618 PRINT#2,"MA";X+GX;",":sY+tGY;";":GOTO 4628 
4619 PRINT#2,"MA"; (X+GX) *MX;","3 (Y+GY) *MX 

462@ IF TL>@® THEN 4678 

4621 IF CC/1@0=INT(CC/18@) THEN 4625 

4622 IF (CC-5)/10=INT((CC-S)/16) THEN 4630 

4623 IF (CC-2.5)/S=INT((CC-2.5)/5) THEN 4635 
4624 GOTO 4670 

4625 ON HP GOTO 4626,4626, 4628, 4629 

4626 PRINT#2,"S1@.18,0.22;PU;PR@,1;":GOTO 4640 
4628 PRINT#2,"MR -1,1;":GOTO 4642 

4629 PRINT#2,"MR"; (-1) *MX;7,";MX:GOTO 4648 

463@ ON HP GOTO 4631,4631,4633,4634 

4631 GOTO 4672@ 

4633 PRINT#2,"MR -1,3;":GOTO 4648 

4634 PRINT#2,"MR"; (-1) *MX;",";3*MX:GOTO 4640 
4635 ON HP GOTO 4636,4636, 4638, 4639 

4636 GOTO 4672 

4638 PRINT#2,"MR -1,5;":GOTO 46490 

4639 PRINT#2,"MR"3 (-1) *MX;",";5*MX 

4648 ON HP GOTO 4641,4641,4643,4644 

4641 PRINT#2,"LB";CC:PRINT#2,CHRS(3):GOTO 4652 
4643 PRINT#2,"LS1;":PRINT#2,"°PL":PRINT#2, CC: PRINT#2,CHRS(13) CHRS(12 
) CHR$(3):GOTO 465@ 

4644 PRINT#2,"SI25,20":PRINT#2,"LA";CC 

4658 TL=TL+1 

4660 GOTO 4615 

4678 2T = 2T + 1 

4680 NEXT Yl 

46980 2T = @ 

470@ FOR Xl = 6 TO 1 STEP .25 

47168 DS = C + D * Xl 

4728 IF DS = @ GOTO 4838 

4730 Yl = (CC - A - B * Xl) / DS 

4748 IF Yl < @ OR Yl > 1 GOTO 4830 

4750 X = X88 + MG* Xl 

4760 Y = YO + MG* Yl 

4770 IF X > PX OR Y > PY GCTO 4828 

4780 IF 2T = @ THEN 4816 

4798 ON HP GOTO 4791,4791,47393,4794 

4791 PRINT#2,"PD;PA";X+GxX;",";Y+GY;";":GOTO 4820 
4793 PRINT#2,"DA";X+GxX;","3%+GY3":":GOTO 48220 
4794 PRINT#2,"DA"; (X+GX) *MX37,"3 (Y+GY) *MX 

480@ GOTO 4828 

4810 ON HP GOTO 4811,4811,4&13,4814 

4811 PRINT#2,"PU;PA";X+Gx;7,";Y+GY;";":GOTO 4828 
4813 PRINT#2,"MA";X+GX;",";Y+Git3";":GOTO 48220 
4814 PRINT#2,"MA"; (X+GX) *MX;","3 (Y+GY) *MX 

4828 ZT = ZT + 1 

4830 NEXT Xl 

4840 NEXT J 

4850 NEXT I 

4860 TL=2 

4870 NEXT CC 

4875 ON HP GOTO 4886,4888,4225,489@ 

4880 PRINT#2,"SP2;S10.2,08.3;":PRINT#2,"PU;PA 16,-3;” 
4881 PRINT#2,"LB":PRINT#2,"FC-Contour Data File - ";FS$;" - Date : "; 
DATES 

4882 PRINT#2, CHRS$(3) 

4883 GOTO 4900 

4885 PRINT#2,"PS2;":PRINT#2,°LS2;":PRINT#2,"MA 18,-3;" 
4886 PRINT#2,"PL":PRINT#2,"°FC-Contour Data File - ";F$;" - Date : "; 
DATES 

4887 PRINT#2,CHRS(13) CHRS’:i2) CHRS(3) 
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4888 GOTO 4900 

4890 PRINT#2,"°SP2":PRINT#2,"S1I40, 30": PRINT#2, "MA"; 10*MX;","73 (-3) *MX 
4891 PRINT#2,"LAPC-Contour Data File —- ";F$;" - Date : ";DATES 
4980 CLS:GOSUB 9800 

4918 LOCATE 10,1@:PRINT “INPUT ‘MORE’ for" 
4928 LOCATE 12,1@0:PRINT "CONTOURS added," 
4930 LOCATE 14,1@:PRINT "*NEWS for fresh plot,” 
4948 LOCATE 16,10:PRINT " 

495@ LOCATE 16,1@0:INPUT "or ‘END'"; MS 

4960 IF M$ = "MORE" GOTO 4230 

4965 CLOSE#2 

4970 IF M$ = "NEW" GOTO 49800 

4980 MS$=" " 

4998 CLS:GOTO 2000 

4995 REM ***** Hewlett-Packard (HP-GL) ****% 
5100 PRINT#2,"PS3;”" 

5105 PRINT#2,"1P170,602,15370,10602;" 

5110 PRINT#2,"SC -5,192,-5,135;" 

5120 PRINT#2,"SP1;" 

5140 PRINT#2,"PU;PA";GX;",";GY;";" 

5141 PRINT#2,"PD;PA";PX;",";GY;";" 

5142 PRINT#2,"PA";PX;",";PY;";" 

5143 PRINT#2,"PA";GxX;",";PY;";" 

5144 PRINT#2,"PA";:GX;",";GY;";" 

5158 FOR I=108 TO (NC-10) STEP 10 

5155 FOR J=1@ TO (NR-10) STEP 10 

5168 PGXO2=GX+(MG*I) 

S161 PGX1=PGX6-2: PGX2=PGX0+2 

5162 PGY8=GY+(MG*J) 

5163 PGY1=PGY@-2: PGY2=PGY@+2 

5170 PRINT#2,"PU;PA";PGX1;",";PGY0;";":PRINT#2,"PD; PA"; PGX2;",";PGYO; 


5175 PRINT#2,"PU;PA";PGXO;",";PGY1;";":PRINT#2,"PD; PA"; PGX@;",";PGY2; 


5188 NEXT J 

5185 NEXT I 

5198 GOTO 4230 

8995 REM ***** pC-Contour Heading ***** 
9860@ PS$="PC-Contour : “":TS=PS+FS$ 

9210 LOCATE 1,2:PRINT TS 

9928 LOCATE 1,30*SR%:PRINT DATES 

9030 PRINT °° s2222E222seee2Se e255 SSS e SESS SSE55See==" 
9848 PRINT 

9850 RETURN 

9500: CLS 

9510 GOSUB 92808 

9528 PRINT: PRINT " Program complete" 
953@ PRINT:PRINT 

9540 END 


3.5 PC-SURFPLOT 


As mentioned at the start of this chapter, a plan or contour drawing presents 
information in 2-dimensions. There are often occasions when there are 
advantages in surface modelling, that is, using an isometric projection to 
display the data. This program can use the same data files, either prepared 
by PC-Grid or PC-Random, and previously used by PC-Contour, to display 
four different views of the area at 90° to one another. In all cases there will be 
some vertical exaggeration to enhance the presentation, and this is 
automatically adjusted to present the display in the centre of the screen. The 
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user has the option of zooming in to enlarge the grid details if necessary. As 
before, screen dumps and plots are available. 


3.5.1 PC-Surfplot — subroutine index 


Line 

numbers Function 
(a) 10-90 Initialization and control 
(b) 100-280 Data file selection 
(c) 1500-1700 Data file read from disc 
(d) 2000-2640 Data file information and screen selection 
(e) 40004330 Graphics screen display 
(f) 5000-5390 Selection of zoom or rotate for graphics 


(g) 6000-6420 Data preparation for plotting 
(h) 7000-8500 Plotter routines 

(i) 9000-9110 Clear screen routine 

(j) 9500-9560 Termination of program 


(a) Initialization and control Line numbers 10-90 


The program Is initialized and the installation file checked for the availability 
of graphics screens and colour for subsequent display of the isometric 
projections. This program requires a Color/Graphics Adapter to be able to 
zoom and rotate the displays. 


(b) Data file selection Line numbers 100-280 


This program gives the user an opportunity to display either Glapwell or 
Glenmore, stored earlier from PC-Grid and PC-Random, or their own data 
file. 


(c) Data file read from disc Line numbers 1500—1700 


The user's data file is entered when requested, the data read from the disc 
and the level of the grid points stored in ZD(NR,NC). 


(d) Data file information and screen selection Line numbers 2000—2640 


The data file information display is similar to that of the previous programs 
in this chapter. The user is then given a choice of displaying the isometric 
projection in colour (medium-resolution) or monochrome (high-resolution) 
on the graphics screen. This program requires a Color/Graphics Adapter 
card to enable the display to be rotated, etc., before plotting. 


(e) Graphics screen display Line numbers 4000-4330 


This display uses the alternative screen presentation of Fig. 1.5, allowing a 
wider plotting area, with two lines of text above and below the graphics 
window. 
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Fig. 3.7. PC-Surfplot screen dump of isometric projection (View 1) of Glenmore on 
Screen 1 (medium-resolution). 


The program presented in this section takes each block of grid points in 
turn and calculates their levels YI-Y4 with X-coordinates X1-X4. and this 
grid square is plotted on the screen. The routine shown uses a hidden line 
algorithm so that old plotting underneath the new grid square is removed 
using the PAINT command. On completion of the mesh diagram. a block 
base 1s drawn showing the vertical height differences at the three visible 
corners of the plot, Fig. 3.7. 


(f) Selection of zoom or rotate for graphics Line numbers 5000-5390 


The first plot of any new data file will be View 1, i.e. a standard 0° isometric 
projection, Fig. 3.8. The user is now given a choice of zooming in or out, or 
of rotating the surface model, Figs 3.9, 3.10 and 3.11. In the zoom mode, 
changes are set to 20% of the original, but this factor can be changed by the 
user if required. 

If the screen presentation Is satisfactory and e for END selected, a further 
list of plotting options is presented: 


Do you wish to save Screen to | 
(1) Disc (2) Printer (3) Plotter (4) End 


thus allowing a screen version to be saved as a picture to disc for future 
recall. us a screen dump to a printer. or as a hardcopy on a plotter. 
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fi ce 
ee ae 

> ‘a 1% 
fade ret 





Fig. 3.8 PC-Surfplot screen dump of isometric projection (View 1) of Glapwell on 
Screen 2 (high-resolution). 








Fig. 3.9 PC-Surfplot screen dump of isometric projection (View 2) of Glapwell on 
Screen 2 thigh-resolution). 
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Fig. 3.10 PC-Surfplot screen dump of isometric projection (View 3) of Glapwell on 
Screen 2 (high-resolution). 
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Fig. 3.11 PC-Surfplot screen dump of isometric projection (View 4) of Glapwell on 
Screen 2 (high-resolution). 
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(g) Data preparation for plotting Line numbers 6000-6420 


The data needs to be rotated before plotting and this routine prepares the 
data depending on the value of ROT%. This routine is used by both the 
screen display and the plotter. 


(h) Plotter routines Line numbers 7000—8500 


Plotter routines are provided for the three plotters discussed earlier, Apple 
410, Epson HI-80 and Hewlett-Packard HP-7475A. The frame is drawn by 
pen 1 (black), the grid mesh by pen 3 (green), the block base by pen 4 (blue) 
and the plot titled by pen 2 (red), as shown in Fig. 3.12. 


(i) Clear screen routine Line numbers 9000-9110 


This routine is used to clear the graphics screen and identify the program and 
isometric projection viewed. 


(j) Termination of program Line numbers 9500-9560 


This is the final display screen and returns the user to BASIC command 
level. 





Fig.3.12 PC-Surfplot isometric projection (View 1) of Glapwell spoil-heap on plotter 
for comparison with Fig. 3.8. 
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3.5.2 PC-Surfplot — numeric variables 





B 

BBY 1-3 
BTY1-3 
BX1-3 
BY 1-3 
DF % 
DISP% 
GCOL 
GROW 
GX 

GY 


MIN 
MX 
NC 
NCOL 
NR 
NROW 
PRTX 
PRTY 
PX 

PY 
ROT % 
S 
SCD&% 
SCR& 
SEL& 
SR& 
SX 

SY 
VD% 
VEX 
X-X4 
XCEN 
Yu-Y4 
YCEN 
YSMIN 
ZH 
Z(58,58) 


ZD(59,58) 


= Number of pixels in isometric grid (horizontal) 


= Variables used to define block base 
= Variables used to define block base 
= Vaniables used to define block base 
= Variables used to define block base 
= Data file section 

= Selected screen save location 

= Graphics screen column variable 

= Graphics screen row variable 


= X-coordinate at lower left-hand corner of plot 
= Y-coordinate at lower left-hand corner of plot 


>Integer variables used in loops 





= Maximum Z-value in data 
= Minimum Z-value in data 


= Multiplication factor used by Epson HI-80 plotter 


= Number of columns in grid 

= Column variable in isometric routine 
= Number of rows in grid 

= Row variable in isometric routine 

= X-location of title on plot 

= Y-location of title on plot 


= X-coordinate at upper right-hand corner of plot 
= Y-coordinate at upper right-hand corner of plot 
== Integer variable to select angle of rotation 

=: Number of pixels in isometric grid (vertical) 


= Integer variable to select screen display 
= Integer variable for screen selection 
= Integer variable to select program 


= Integer variable for graphics screen resolution 


== Horizontal display width in user units 
== Vertical display height in user units 
= Video display variable 

= Vertical exaggeration 

= .-coordinates of grid points 

= X-coordinate of centre of plot 

= Y-coordinate of grid points 

= Y-coordinate of centre of plot 

== Minimum value of Z on screen 

= Read Z value of grid points 

= Array to store Z-values of grid points 
= Array to store Z-values of grid points 
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3.5.3 PC-Surfplot — string variables 





A$ = INKEY$ 

DATES Vaniable used to return the date of the computer’s calendar/clock. 
in US format MM-DD-YYYY 

F$ File name for data file 

Q$ Question response (Y/N) or (y/n) 

SCR$ File name to save picture of screen on disc 

SCRPIC$ = SCR$ + “.PIC” 

TIMES Variable used to return the time of the built-in clock in 
HH:MM:SSS format 

V$ Variable to store view shown in isometric projection 





3.5.4 PC-Surfplot — BASIC program 


10 
12 
14 
16 
20 
Col 
22 
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sc 
220 


REM <SURFPLOT> Program for IBM PC and Compatibles 
REM Version 81-81-1987 
REM (C) Copyright 1987 P.H.Milne 
REM All rights reserved 
REM Output available to <1l> Apple 410, <2> Epson HI-88, <3> HPGL 
our Plotter 
REM RS-232C Setting in Line 4400 - OPEN "COM1:9600,E,7,1,CS2000,DS 
@@" FOR OUTPUT AS #2 
REM Can access data files on disc ZOOM and ROTATE display 
REM This version includes Graphics printer dump of Screen 
CLEAR 
PCFS="PCSCREEN" 
OPEN *1",4#3,2CFS 
INPUT#3,SCR% 
INPUT#3,SCD% 
CLOSE#3 
ON (SCR%+1) GOTO 62,64,66 
SCREEN @:WIDTH 48:GOTO &@ 
SCREEN 1:GOTO 68 
SCREEN 1:GOTO 86 
IF SCD%=1 THEN 880 
COLOR 8,8 
KEY OFF 
DIM 2(59,50), 2D(58,58) 
SR%=1:DF%=0:DISP%=1: ROT%=1:VD%=O8:FS=" "* 
CLS:GCSUB 9800 
PRINT " Seesesssessesesessessessssseeseeseesse="3spRINT 
PRINT " THIS PROGRAM REACS INFORMATION FROM " 
PRINT " REGULAR SQUARE GRID POINTS INTO A ®” 
PRINT " DATA FILE FOR PLOTTING BY PC-SurfPlot" 
PRINT 
PRINT " To Continue, Press <Return> ":WHILE RS="":RS=INKEYS :WEND 
PRINT: PRINT 
PRINT " Select Data File for Display" 
PRINT: PRINT " <1> Glapwell, <2> Glenmore or":PRINT " <3> Other Gi 
file";: INPUT DF% 
IF DF%<1 OR DF&>3 THEN 208 
PRINT: PRINT 
ON DF% GOTO 25@,2708,1588 
FS="Glapwell" 
GOTO 15280 
FS="Glenmore" 
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2808 GOTO 1520 

1495 REM ***** Data File Disc Routine ***e# 

1508 PRINT:PRINT " Enter DATA File Name "; 

151@ INPUT FS 

1526 PRINT:PRINT " One moment please .... reading data" 
1530 OPEN “I",#@1,FS 

1540 INPUT#l, NR 

1558 INPUT#@1, NC 

1560 INPUT#1, MAX 

1578 INPUT#1, MIN 

158@ INPUT#1, DX 

1598 INPUT#1, DY 

1600 FOR I#9@ TO (NR-1]) 

1610 FOR J#8 TO (NC-1) 

1620 INPUT#1, ZH 

1638 Z2ZD(I,J)#ZH 

1648 NEXT J 

1650 NEXT I 

16608 CLOSE@] 

1678 PRINT 

1680 PRINT *® Do you wish to check Data File (y/n)";: INPUT Q$ 
16980 IF Q$#"n" THEN 2120 

17800 IF QS<>"y" THEN 1670 

1995 REM ***** Data File Information ***** 

2000 CLS:SCREEN 1 

20108 GOSUB 9808 

2020 PRINT  sesstee2rsee2eee2es225S 222222222 222222===" 
2030 PRINT: PRINT: PRINT " THE DATA FILE - “";FS:PRINT: PRINT ™ HAS NOW B 
EEN READ, GIVING :-" 

2040 PRINT 


";DY/(NR-1);" Units" 
"sDX/(NC-1);" Units" 


2050 PRINT ™" No. Rows = ";NR;" Now Cols. = “;NC 
2060 PRINT " Minimum value = ";MIN 

2070 PRINT " Maximum value = ";MAX 

2080 PRINT " Range = "; (MAX-MIN) 


2090 PRINT 1 Map row 
21008 PRINT “ 1 Map column 
2110 IF QS="y" THEN 2148 
2128 CLS:GOSUB 9000 

2130 PRINT " SSSsssssssessessseessscssssesse2=s=====" 

2140 PRINT: PRINT 

215@ PRINT ™ Do you wish to";TAB(18);"<l> Plot Surface":FRINT TAB(L2) 
»*<2> Run PC=Contouc * 

2160 PRIGT TAB(1a);"<3> Quit *;:sINPUT SEL 

2i70 IF SELa<l OR SELG>3 THEN 2146 

2180 If SEts=3 THES 9509 

2190 IF SEL%=1 THEN 2508 

2298 LCAD “A:CCNTOUR",R 

2495 REM *e*** Screen Selection. ****s 

2500 CLS 

2510 GOS UE 9200 

2520 PRINT " sewer tesetetsteteteeeeseesere2eee=2=22==== 
2530 PRIST: PRINT THS(3);"Select SCREEN";TAB(20);"<1> Colour -" 


25408 PRINT TAB(23);" Low resolution" 
2550 PRIST:PRINT TAS(20);"<2> Slack & White -" 
2560 PRINT TAB(2¢2);" Hian resolution" 


2570 PRINT TAE( 23) 92 1HPUT SRe 

2580 IF SAt<) OR SRI>2 THEN 2530 

2596 IF RCT%=]1 THEN NCOL=NC-1] :NROW=NR-1 

2600 PRINT: PRINT: PRINT “ One moment please .... preparing cata" 
2610 IF SF&%=1 THEN 6000 

2620 SCKEEN 2 

2630 YCEN=aN:XCEN=8:SX=160*SRB&:S7Y=100 

2640 GOTO 6080 

3995 REM *S**" Screen Plotting Routine ***** 

a@O0 DISC e2l: YCENsd: XCEN=0:°%=2=160*SR3:SY=lC0 


seca mrpetys yy "Sc SAR Qe 
twas Lk e : THES, rs es, Pues "“Gos%=3 


4070 
40880 
4090 
4100 
4116 
4120 
4138 
4140 
415@ 
4160 
4178 
4180 
4190 
4200 
4218 
4220 
4230 
4248 
4258 
4260 
4270 
4286 


4290 
4360 
1310 


4320 
4330 
+395 
4400 
4405 
+410 
4420 
4430 
4440 
4450 
1466 
4470 
3490 
4499 
4566 
4582 
4505 
1510 
“S15 
4916 
4517 
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VIEW SCREEN (8,18)-(310*SR%,180) ,,RGBSE 
WINDOW (XCEN-SX,YCEN=-SY) -(XCEN+SX,YCEN+SY) 
VEX=5@/ (MAX-MIN) : YSMIN=INT( (MIN*®* VEX) /1@) *10+5 
X8=-10*SRB%: YG=(70-YSMIN) 

CLS 

GOSUB 90888 

S=2.25:B=3*SR8% 

IF ROT%=2 OR ROT%=4 THEN XG=10*SR% 

FOR GCOL=1 TO NCOL 

FOR GROW=1 TO NROW 

X1=X0+ (GCOL-GROW) *B 

Y1l=Y0@+Z (GROW,GCOL) *VEX-(GROW+GCOL) *S 
X2=X1-B 

Y2=Y@0+Z (GROW ,GCOL=-1) *VEX- (GROW+GCOL=1) *S 
X3=X1+B 

Y3=Y0+Z (GROW-1,GCOL) *VEX=(GROW-1+GCOL) *S 
X4=X1 

Y4=Y0+Z (GROW-1 ,GCOL=1) *VEX-(GROW=-1+GCOL=-1) *S 
LINE (X1,Y1)-(X2,Y2),2 

LINE -(X4,Y4),2:LINE -(X3,Y3),2 

LINE -(X1,Y1),2 

Y=(Y1+Y4)/2 

PAINT (X4,Y),0,2 

LINE (X1,Y1)-(X2,Y2) ,2 

LINE —-(X4,Y4),1:LINE -(X3,Y3),1 

LINE -(X1,Y1),1 

NEXT GROW,GCOL 


123 


LINE (X0+(-NRCW) *8,YO-(NRCW) *S+Z (NRCW,@) *VEX) -(X0+(-NROW) *B,YO-( 
NROW) *S+YSMIN) 
LINE (XO0+(NCOL-NRCW) *8, r%0-(NCOL+NROW) *S+Z (NROW,NCOL) * VEX) -(XO+(N 
COL-NROW) *B, Y@-(NCOL+NRCW) *9+YSMIN) 
LINE (X0+(NCOL) *8,YC=(NCOL) *S+Z(0,NCCL) * VEX) = (X04+(NCOL) *8,YO=(NC 
JL) *StyYonu iN) 
CINE (Xv+(-NROW) *8,YC-{NRCwW) *S+YSMIN) -(X6+(NCOL-NROW) *8,YC-(NCOL 
+NROW) *S+YSMIN) 


LINE -—(A0+(NCCL)*8,Y¥O=—45COL) *S+YSMIS) 


GOTO Svvo 

REM ***** Plotter Routine *x*r* 

OPEN "COM1:960U0,2£,7,1,c22008,0Si50060" FOR OUTEUT AS 32 
, 4 

CIS P%t=2 


IF SR*&=2 THEN Sk=SX/2 

IF SR*€=2 THEN B=6/2 

Ir ROT%=2 OR ROTS=4 THES Kd=10 

VEX=56/ (MAX<MIN) :YSMINSIST ( (MIN* VEX) /12) *1le+5 
GK=XCEN=-SX :GY=YCEN=-SY: EX=KCEN+SX3 PY=YCEN+SY sMX=2670/ 320 
FRTX=GX+2: PRTY=GY+2 

If ROT%=1 THEN COLENC-L:hRCWH=NR-] 

OM HP® GOTO 45¢¢,4580,4520,4575 

REM rarrkne HP=l, (EPGCL} Ferran 
PRINT#2,"PS3;" 

PRINT# 2, "IP 170,682,1S272,)0602:" 
PRINT#2,"SC"°GX-" ,“ sere", "eGo" We pye" 2 
PRINT#2,"SP1;" 
PRINT#2, "PU; PA" SSK 7","232e";" 
PRINTS2,"P0O; PA": 2x3", "5Ge3"5™ 

PRINT# 2,5 PA SPR", 2 Peete 

PRINT# 2," PA**GXe ny" eres” 

PRINT# 2; PA" 2GxXe"," 2Gerts” 

PRINT#2Z,; S33" 

GOTO 6Ue¥u 

REM ***** Hp=3, Acple gla *#ees 
PRINT#2,°SP33" 

PRINT# 2y WO" 7OXe  p PGE ep PPR Eg PPR es” 
PRIATH 2; PSa<" 
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4565 PRINT#2,"MA";GX;",";GY;";":sPRINTH2,"DA";PX;","3GY;";":PRINT#2,"D 
Ane PX ee FT PYe "2 "2 PRINT#H2, "DA" ;GX;","sP¥2";"2PRINT#H#2, "DA" 7GXi","2GYi": 


4570 PRINT#2,"PS3;":GOTO 60080 

4574 REM ***** Hp=4, Epson ***#* 

4575 PRINT#2,"IWO,8,2670,1920" 

4580 PRINT#2,"MA1335,960":PRINT#2,"OR" 
4581 PRINT#2,"SPl1":PRINT#2, "LTO" 

4582 PRINT#2,"MA" ;GX*®*MX;"," :GY*MX 

4583 PRINT#2,"DA";PX*®*MX;",";GY*MX 
4584 PRINT#2,"DA";PX*MX;",";PY*MX 
4585 PRINT#2,"DA" ;GX*MX;","3 PY®MX 
4586 PRINT#2,"DA";GX*MX;","sGY*MX 
4598 PRINT#2,"SP3" 

4595 REM ***** Plotter Output ****s 

4600 FOR GCOL 6 TO NCOL 

4610 FOR GROW 6 TO (NROW-1) 

4620 X1=X0+(GCOL-GROW) *B 

4625 Yl=Y@+Z(GROW,GCOL) *VEX-—(GROW+GCOL) *S 

4630 X2=Xl1-B:Y2=YO+Z (GROW+1,GCOL) *VEX-(GROW+1+GCOL) *S 
4640 IF GROW<>@ THEN 4660 

4650 ON HP® GOTO 4651,4651,4653,4654 

4651 PRINT&2,"PU;PA";X1;3","7Y1;5";"%:GOTO 4668 

4653 PRINT#2,"MA";X1:",";Y1;"3":GOTO 4660 

4654 PRINT#2,"MA";Xl1*MX;3"," 5; Y1*MX 

4660 ON HP% GOTO 4661,4661,4663,4664 

4661 PRINT#2,"PD;PA";X2;",";Y2;";":GOTO 4670 

4663 PRINT#H2,"DA";X2;","sY¥23";":GOTO 4678 

4664 PRINT#2,"DA";X2*MX3",";Y2*MX 

4670 NEXT GROW,GCOL 

4680 EFCR GROW = JD TO NROW 

4698 FOR GCOL = J TO (NCOL-1) 

A700 Xl1=xXC+(GCOL-GROW) *B 

4710 YlaYo+Z7Z (GRCw,GCOL) *VEX=(GRCwW+GCCL) *S 

$720 X2=Xcd+(GCCL+1-GRCw) *B 

4730 Y2=2"%0+Z (GROw,GCOL+]) *VEX-(GROW+GCOL+1) *S 

4740 IF GCOL<>¢ FHEN 4768 

4753) @M FPt. GOTO 4$751,4751,4/53,48 754 

475] PREP? , “PS Pan exe ele se * GCs 4760 

4753 PRs e ee he Seles eV 1S se 2GCio <6 

R7S4 ERINTEZ MA" XL", VY ISMX 

4768 ON =P3 GOTO 4761,4761,4763,4764 

676) PRENTSE2, "Fes SA" eKo "eV 2s" es" GCS 4770 

C76 ESL Nee 2, “Sr Ko eZee GOTO 477 
4764 PREST32,"fa":K2*MX2"," 2 Y2*MX 

4770 NEXT GCOL,GRCW 

4775 BAL=AG+ (-MACH) *BS ETYL=AYO-(NRCW) *f*+7 (NROW,C) *VEX: BBYL=YO-(NROW) *S 
+YSMIN 

4786 BL Z2=Ko+ (NCTCL-NROW) *BsBTY2=YO-(NCOL+NROW) *S+Z (NROW ,NCGL) *VEX: 5572 
=YO0-(NCTL ARON) FSETYSMIN 

4785 BAlzL e+ (CCL) *BsBTY3=YO-(NCOL) *$+2(0,NCOL) *VEX: BBY3=YC=—(NCOL) *S+ 
YSMIN 

4790 O28 Hrs GOTT 4800,480U0,4660,492¢ 

4799 S**8* Host, LHP-GL) *****# 

4600 PRENTE2,"“Sr4:" 

4801 PRINT32,"FU;PA";8Xi:",";BTYis"s" 

“oo. Pectacy rerraA goes “2° Sey 1s *s” 

4810 PRInT#2,"PA";3X2:","sEBY2; "3" 

4815 PRInT$2,"PA";8X3;","s;BBY3;"3" 

4820 PRINT#2,"PhA";BX3;",";BTY3;";PU;" 

4825 PRINT#2,"PA";3X2;",";BTY2;"5" 

4830 PRE:T#2, "Pl; FA";BX2;",";BBY2;";PU;" 

£035 Pears 2;°SPz2i" 2 PRINTH2, "S10.2,0.22" 

484€ PRINTA2,"FU;FPA ";PRTX;",";PRTY;";" 

4845 PRIn7T32,"Le";"PC-SurfPlot - "Skor" = View "sVS." = Date » “ZEA TE 


a) 
u 
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4850 FRINT#2, CHRS(3) 
4855 GOTO 4989 

4859 REM ***** HP=3, Apple 410 ***#*# 
4860 PRINT#2,"MA";BX1;",";BTY1; 


4865 PRINT#2,"DA";BX1; “ie co ee 
4870 PRINT#2,"DA"; BX23;","; BBY2> "3" 
4875 PRINT#2,"DA";BX3;",";BBY3;";" 
4880 PRINT#2,°CA" ;BX3;","sBTY33";" 
4885 PRINT#2,"MA";BX2; mG oes 


4890 PRINT#2,"DA";BX2; 7BBY2: 

4895 PRINT#2,"PS2:":PRINT#2,"LS3; 

4900 Se a ee 

4905 PRINT#2,"PL":PRINT#2,"PC-SurfPlot -— ";FS$;" 
Ss DATES 

4910 PRINT#2, CHRS$(13) CHRS$(10) CHRS(3) 

4915 GOTO 49898 

4919 REM ***** Hp=4, Epson ***** 

4920 PRINT&#2,"MA" ;BX1*MX;",";BTY1]*MX 

4925 PRINT#2,"CA";BX1*MX;",";BBY1*MX 

4930 PRINT#2,"DA";BX2*MX;",";BBY2*MX 

4935 PRINT#2,"DA";BX3*MX;",";BBY3*MX 

4940 PRINT#2,"DA";BX3*MX;",";BTY3*MX 

4945 PRINT#2,"MA";BX2*MX;",";BTY2*MX 

4950 PRINT#2,"CA";BX2*MX;",";BBY2*MX 

4955 PRINT#2,"SP2":PRINT#2,"S1460, 30" 

4960 PRINT#2,"MA": (GX®*MX4+10)3",73 (GY*MX+5) 

4965 PRINT#2,"LA":PRINT#2,"PC-SurfPlot ~ ";FS$;" 
“sDATES 

$980 CLOSE#2 

4990 IF SR%=2 THEN SX=SX*2:B=B*2 

4995 REM ****#® Selection of ZOOM or ROTATE ***** 


S5u0@ ESEP:DISP$=1: LOCATE 24,2:PRINT “ZOOM = z, ROTATE = cr, 


SU10 ASHENKEYS 
$020 IF *S=""_ TREN GOTO 501d 


SU30 IF iIeSTR("zre*,AS)=0 THEN GOTO Stee 

SJ4¢ IF i¢="2" THEN $136 

USO iF 2S="e" THEN 527¢ 
eon Ile es ae "CCTATE = Select anale cz 
aR PUREE <i> B ,<2>.90,<¢2> 180,<4> 270,<5> ERE"; 
OUT AS=ERKEYS: IF As="" TEEN $876 


2O@C £F INSTR("12345",AS)=0 THEN S@7e 
Se9C RCTA?SVAL(AS) 
Saul LOCREE 24,22PRINT “ 


SL1lg if 2OTt=5. THEN 9500 
Sl20 COPD. 6e09 


$13@ LOCATE 24,2:PRINT “ZOGM IN=i, OUT=o, 
S,2:PRINT"uUPsu, DCWN=d, RiGHT=r, LEFT=1, END=e"; 
$143 AS=INKEYS 

S150 IF As="" THEN 5149 

SIG? IF SSSTR(*ioudelr",AS)=0 THEN GOTO 513¢ 
S170 LOCATE 24,2:PRINT "™ 

Speirs Se 

Sige rF +S="e" THEN 5272 

5i98 £F +5="0" THEN SX=S¥/.2:SY=SY/.8 

S2vG IF 2+S="i" THEN SX=SX*.,é:SY=SY*.8 

9226 IF AS="r" THEN XCEN=XCEN-SX*, 2 

922@ IF AS="1" THEN XCEN=XCEN+SX®*. 2 

9238 IF AS="u" THEN YCEN=YCEN-SY*, 2 

5248 IF 2$="d" THEN YCEN=YCEN+SY®*, 2 

5253 wintGw (XCEN-SX,YCEN-SY) -(XCEN+SX, YCEN+SY) 
2268 GGTS 4382 

227g UCLare. 24; 22PRINT 


t 


- View 


- View 


rotation 


"sVS7" 


-V53” 


~ Date 


NO 


126 SURFACE MODELLING 


5288 LOCATE 24,2:PRINT "Do you wish to save Screen to";:LOCATE 25,2:P 
RINT"<1>Disc<2>Plotter<3>Printer<4>End"; 
5298 INPUT PLOT@:IF PLOT&<1 OR PLOT8>4 THEN 5278 


§380 LOCATE 24,2:PRINT ® "3: LOCATE 
25,2: PRINT" "ys 

5310 ON PLOT% GOTO 5326,5488,5588,9508 

§328 LOCATE 24,2:PRINT "* ";:LOCATE 2 
5,2: PRINT" "4 


5330 LOCATE 24,2:PRINT "Enter File Name for saved poreen ";:LOCATE 25 
,2: PRINT"? (8 letters max.)"*; 

5348 INPUT SCRS:IF SCRS$="" THEN 53290 

5358 DEF SEG = &HB886 

5368 SCRPIC$#SCR$+". PIC® 

5378 BSAVE SCRPICS,8,&H4660 


53868 LOCATE 24,2:PRINT ” ";:LOCATE 2 
5,2: PRINT® "" 

5398 GOTO 5868 

5488 LOCATE 24,2:PRINT * , “;: LOCATE 
25,2: PRINT® -% 


5410 LOCATE 24, 2:PRINT "Select Plotter -<l> Apple, <> Epson"; : LOCATE 
25,2:PRINT “or <3> Hewlett-Packard (HP-GL) "; 

5420 INPUT PX%:IF PX%&<1 OR PX%>3 THEN 5400 

54380 LOCATE 24,2:PRINT * ";: LOCATE 
25,2: PRINT” : 
5440 IF PX*=2=]1 THEN HP%=3 

5458 IF PX*22 THEN HP%=4 

5460 IF PX%=3 THEN HP%=] 

5478 GOTO 4488 

5495 REM ***** Graphics Dump *****# 

5588 LCOPY 

5518 GOTO 5602 

5995 REM ***** Data preparation for Plotting ***#* 

6008 VD%t=1 

6810 ON ROT% GOTO 6028,61908,6210,6338,9500 

6828 NCOL=NC-1:NROW=NR-)] 

6830 FOR I = 8 TO NROW 

6848 FOR J = @ TO NCOL 

60560 2Z2(I,J)=ZD(I,J) 

6060 NEXT J 

60708 NEXT I 

608@ IF DISP%=1 THEN 400@ 

6098 IF DISP%=2 THEN 468690 

6188 NCOL=NR-1:NROW=NC-]1 

6116 K=NCOL 


=e 


6128 FOR I = @ TO NROW 
6138 FOR J = @ TO NCOL 
6140 Z(1I,J)=ZD(K,1) 
615@ K=K-1 


6160 NEXT J 

61706 K=NCOL 

6188 NEXT I 

6198 IF DISP%=1 THEN 48586 
6200 IF DISP%=2 THEN 4680 
6218 NCOL=NC-1:NROW=NR-1 
6228 K=NROW: L=NCOL 

6230 FOR I=@ TO NROW 

6248 FOR J=@ TO NCOL 

6258 2(1,J) =ZD(K,L) 

6268 L=L-1 

6278 NEXT J 

6280 K=K-1 

6298 L=NCOL 

6300 NEXT I 

6310 IF DISP%=1 THEN 4956 
632@ IF DISP%=2 THEN 4600 
6338 NCOL=NR-1:NROW=NC-] 


6348 
6358 
6368 
6370 
6388 
6398 
6408 
6418 
6426 
8995 
9008 
9916 
98620 
9636 
9640 
9858 
9668 
9870 
9880 
9698 
9180 
91190 
9495 
9500 
9518 
9528 
9538 
9548 


PC-SURFPLOT 


K=NROW 
FOR I = @ TO NROW 
FOR J = @ TO NCOL 
Z(I,J)=ZD(J,K) 
NEXT J 
K=K-1 
NEXT I 
IF DISP%=1]1 THEN 4858 
IF DISP%=2 THEN 46808 
REM ***** pC-SurfPlot Heading ***** 
LOCATE 1,2:PRINT "PC-SurfPlot” 
LOCATE 1,31*SR&:PRINT DATES 
LOCATE 2,2:PRINT FS 
IF VD%=8 THEN 9188 
IF ROT%=1 THEN V$="l: @ Degs" 
IF ROT%=2 THEN V$="2: 98 Degs” 
IF ROT%=3 THEN V$="3:188 Degs" 
IF ROT%=4 THEN V$="4:278 Degs” 
IF ROT%=5 THEN vS$=" " 
LOCATE 2,12*SR&:PRINT “View °;VS 
LOCATE 2,31*SR&*:PRINT TIMES 
RETURN 
REM ***** Program End *****# 
CLS: VD%=9 
LOCATE 2,2:PRINT " 
GOSUB 9890 
VIEW SCREEN (8,18)-(318*SR%,180),,@ 
LOCATE 24,2:PRINT * 


§,2:PRINT” ; 
LOCATE 8,5:PRINT * Program Complete":PRINT:PRINT: PRINT: PRINT 


9558 
9568 


END 
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";:LOCATE 2 


